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ABSTRACT 

We study the physical processes that affect the ahgnment of grains subject to ra- 
diative torques (RATs). To describe the action of RATs, we use the analytical model 
(AMO) of RATs introduced in Paper I. We focus our discussion on the alignment 
by anisotropic radiation flux with respect to magnetic field, which defines the axis of 
grain Larmor precession. Such an alignment does not invoke paramagnetic dissipation 
(i.e. Davis-Greenstein mechanism), but, nevertheless, grains tend to be aligned with 
long axes perpendicular to the magnetic field. When we account for thermal fiuctu- 
ations within grain material, we show that for grains, which are characterized by a 
triaxial ellipsoid of inertia, the zero- J attractor point obtained in our earlier study 
develops into a low-J attractor point. The value of angular momentum at the low-J 
attractor point is the order of the thermal angular momentum corresponding to the 
grain temperature. We show that, for situations when the direction of radiative flux 
was nearly perpendicular to a magnetic field, the alignment of grains with long axes 
parallel to the magnetic field (i.e. "wrong alignment") reported in Paper I, disappears 
in the presence of thermal fluctuations. Thus, all grains are aligned with their long 
axes perpendicular to the magnetic fleld. To gain insight into the origin and stability 
of the low-J attractor points, we study the dynamics of the crossovers in the pres- 
ence of both RATs and thermal fluctuations. We study effects of stochastic gaseous 
bombardment and show that gaseous bombardment can drive grains from low-J to 
high-J attractor points in cases when the high-J attractor points are present. As the 
alignment of grain axes with respect to angular momentum is higher for higher values 
of J, counter-intuitively, gaseous bombardment can increase the degree of grain align- 
ment in respect to the magnetic field. We also study the effects of torques induced 
by H2 formation and show that they can change the value of angular momentum at 
high-J attractor point, but marginally affect the value of angular momentum at low-J 
attractor points. We compare the AMO results with those obtained using the direct 
numerical calculations of RATs acting upon irregular grains and validate the use of 
the AMO for realistic situations of RAT alignment. 
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1 INTRODUCTION 

Polarization from absorption and emission by aligned grains 
is widely believed to trace magnetic field topology (see Hilde- 
brand et al. 2000, Hildebrand 2002; Aitken et al. 2002; 
Hough et al. 1989). The reliability of the interpretation of 
polarization maps in terms of magnetic fields depends cru- 
cially on the understanding of grain alignment theory. 

When, nearly 60 years ago, the polarization of starlight 
was discovered (Hall 1949; Hilner 1949), it was immediately 
explained by absorption by elongated dust grains, which are 
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afigned with respect to the interstellar magnetic field. Since 
then, the problem of grain alignment has been addressed by 
many authors (see a recent review by Lazarian 2007 and ref- 
erences therein). As a result, substantial progress has been 
made towards understanding how these tiny particles can 
become aligned. 

The original theory of grain alignment formulated by 
Davis & Greenstein (1951) is based on the paramagnetic 
dissipation of energy for a grain rotating in an external mag- 
netic field. However, the paramagnetic relaxation time-scale 
for a typical interstellar field strength is long compared to 
the gas damping time. Moreover, this mechanism has been 
found to be more efficient for small grains (see Lazarian 
1997a; Roberge & Lazarian 1999), contrasting with obser- 
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vation data, which have testified that small grains are not 
aligned (Kim & Martin 1995; Andersson & Potter 2007). 

Study of the paramagnetic alignment mechanism was 
given new life by Purcell (1979). In his classical paper, Pur- 
cell suggested three processes that can spin a grain up to 
suprathermal rates (i.e., much higher than thermal value): 
H2 formation, photoemission, and the variation of the ac- 
commodation coefficient over the grain surface. H2 forma- 
tion was identified as the most powerful of these three 
processes. In terms of alignment, fast rotation is advanta- 
geous, as fast rotating grains are immune to randomiza- 
tion by gas bombardment. As a result, paramagnetic dis- 
sipation can give rise to good alignment of angular momen- 
tum with the magnetic field. An obvious limitation of the 
Purcell mechanism is that the spin-up process is most ef- 
ficient, provided that a sufficient fraction of hydrogen is in 
atomic form. Therefore, this mechanism would fail in dark 
molecular clouds where most of the hydrogen is in molecular 
form. In fact, all spin-up processes fail in molecular clouds 
(see Lazarian, Goodman & Myers 1997; Roberge & Lazarian 
1999), while observations have testified that the grains were 
aligned there (see Ward-Thompson et al. 2000). 

Nevertheless, combining the ideas of suprathermal ro- 
tation with magnetic inclusions, in the spirit of the classical 
Jones & Spitzer (1967) paper, researchers hoped to explain 
the observational data. For instance, a model by Mathis 
(1986) provided a good fit to the observed Serkowski po- 
larization curve (Serkowski 1975). Problems concerning the 
Purcell (1979) alignment mechanism became obvious, how- 
ever, when Lazarian & Draine (1999a) reported the effect 
of thermal trapping. This effect stems from the coupling of 
vibrational and rotational degrees of freedom via the inter- 
nal relaxation (e.g. Barnett relaxation; Purcell 1979). Ther- 
mally trapped grains undergo fast flips that average out 
un-compensated Purcell torques. Later, Lazarian & Draine 
(1999b) reported a new effect which they termed nuclear re- 
laxatiorQ; they showed this to be 10^ times stronger than the 
Barnett relaxation. After this, it became clear that most of 
the interstellar grains cannot be spun up by Purcell torques. 

The mechanical alignment mechanism, which is based 
on the relative motion of gas and dust, was pioneered 
by Gold (1952a). However, this mechanism, as well as its 
more sophisticated cousins (Lazarian 1995, 1997; Lazarian 
& Efroimsky 1996; Lazarian, Efroimsky & Ozik 1996), re- 
quires supersonic motion of gas relative to dust (see Purcell 
1969; Lazarian 1994; Roberge, Hanany & Messinger 1995; 
Lazarian 1997b), and this is only available in special cir- 
cumstances. In his original paper, Gold (1952b) proposed 
collisions between clouds as a way to drive supersonic mo- 
tion. Soon after this, Davis (1955) pointed out that such a 
process can only align grains in a tiny fraction of the in- 
terstellar medium (ISM). Other, more promising means of 
obtaining supersonic drift are ambipolar diffusion (Roberge 
& Hanany 1990, Roberge, Hanany & Messinger 1995) and 
MHD turbulence (Lazarian 1994; Lazarian & Yan 2002; Yan 
& Lazarian 2003; Yan, Lazarian & Draine 2004). Neverthe- 



^ The Barnett relaxation arises from the unpaired electron spins 
in the grain. The nuclear relaxation is resulted from unpaired 
nuclear spins. The mechanisms are different and should not be 
mixed up. 



less, while applicable for particular environments, the me- 
chanical mechanisms are unable to explain the ubiquitous 
alignment of dust in a diffuse medium and molecular clouds. 
A more promising mechanism is the mechanical alignment of 
helical grains, first mentioned in Lazarian (1995) and Lazar- 
ian, Goodman & Myers (1997). We briefly described this in 
Lazarian (2007) and Lazarian & Hoang (2007a). The conse- 
quences of this have to be evaluated, but these are beyond 
the scope of the present paper. 

Alignment by radiative torques (RATs) has recently 
become a favored mechanism to explain grain alignment. 
This mechanism was initially proposed by Dolginov & My- 
traphanov (1976), but was not sufficiently appreciated at 
the time of its introduction. Draine & Weingartner (1996, 
1997, hereafter DW96 and DW97) reinvigorated the study 
of the RAT mechanism by showing that RATs induced by 
anisotropic radiation upon rather arbitrarily shaped grains 
can spin-up and directly align them with the magnetic field. 
These papers refocused attention on the RAT mechanism 
and made it a promising candidate to change the grain 
alignment paradigm. It is also very noble that Bruce Draine 
has kindly modified the publicly available DDSCAT code 
(Draine & Flatau 1994) to include RATs. This has enabled 
other researchers to access a useful tool for studying the 
RAT alignment. 

The RAT mechanism seems to be able to address major 
puzzles presented by observations. For instance, the observa- 
tion of polarized emission emanating from starless cores (see 
Ward- Thompson 2000, 2002) initially seemed completely 
unexplainable0 Indeed, all mechanisms seemed to fail in 
such cores, which are presumably close to thermodynamic 
equilibrium (see Lazarian, Goodman & Myers 1997). The 
RATs seem to be too weak as well (DW96). However, Cho 
& Lazarian (2005, hereafter CL05) found that the efficiency 
of RATs increases fast with grain size (see also Lazarian & 
Hoang 2007a), and thus large grains can stiU be aligned in 
dark clouds. They found that grains as large as 0.6/im can 
be aligned in dark clouds by radiation attenuated by the 
column density with Av ~ 10. This was for much higher 
extinction than expected in the absence of embedded stars 
(see DW96), and was claimed to be observed with optical 
and near infrared polarimetry (Acre et al. 1998). 

The pre-stellar cores studied in Ward-Thompson et al. 
(2000) correspond to Ay = 50 — 60; the shielding column, 
assuming uniformity, is approximately half of these values. 
However, Crutcher et al. (2004) pointed out that polariza- 
tion data do not sample the innermost core region^, and 
this provides an explanation for polarization "holes" (see 
Matthews & Wilson 2000, 2002; Matthews et al. 2001; Lai 
et al. 2002). The reported decrease in the percentage po- 
larization with the optical depth also agrees well with the 
findings in CL05. 

The approach of CL05 was further elaborated in the 
studies by Pelkonen et al (2007) and BetheU et al. (2007), 
in which the synthetic maps obtained via MHD simulations 

^ These findings are also contrasting with observational claims 
based on visible and near-infrared radiation (Goodman et al. 
1995). The difference in results was explained in Lazarian (2003) 
The peaks Ay of 150 were claimed for the clouds in Pagani et 
al. (2004). According to Crutcher (2007, private communication) 
these peaks are likely not to produce polarized dust emission. 
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were analyzed. In particular, Bethell et al. (2007) obtained 
the polarization maps were for a simulated cloud with peaks 
of extinction Ay as high as 25. The study showed non-trivial 
nature of radiative transfer in a fractal turbulent cloud, and 
confirmed the the ability of aligned grains to trace a mag- 
netic field in dark clouds and proved a decrease of percentage 
polarization at the highest Ay - Note that the studies above 
were in contrast to the previous studies, which used rather 
arbitrary criteria (e.g. Av = 3) for the alignment to shut 
down, or even more unrealistic assumption that all grains 
were perfectly aligned. 

Similarly, the change of the degree of polarization with 
optical depth observed by Whittet et al. (2001) was also 
explained by RATs in Lazarian (2003, 2007). Detailed mod- 
eling of this effect is provided in Whittet et al. (2007) and 
Hoang & Lazarian (in preparation). 

The encouraging correspondence of the theoretical ex- 
pectationfl and the observed polarization arising from 
aligned grains makes it essential to understand the reason 
why RATs align grains. In our first paper (Lazarian & Hoang 
2007a, hereafter Paper I) , we subjected to scrutiny the prop- 
erties of RATs. Using a simple analytical model (AMO) of 
a helical grain we studied the properties of RATs and the 
alignment driven by such RATs. The analytical results were 
found to be in good correspondence with numerical calcula- 
tions for irregular grains obtained with DDSCAT. Evoking 
the generic properties of the RAT components, we explained 
the RAT alignment of grains in both the absence and pres- 
ence of magnetic fields. Intentionally, for the sake of simplic- 
ity, in Paper I we studied a simplified dynamical model to 
demonstrate the effect of the RAT alignment. This model 
disregards the wobbling of grain axes with respect to the 
angular momentum that arises from thermal fluctuations 
(Lazarian 1994; Lazarian & Roberge 1997), and thermal 
flipping of grains (Lazarian & Draine 1999a, b). The sim- 
plifications allowed us to provide a vivid illustration of the 
RAT alignment, however, there is still a question concerning 
the modification of the alignment by thermal fluctuations, 
as well as the action of additional (e.g. random) torques. 

The first study to combine the RAT alignment and the 
physics of thermal fluctuations and flipping was carried out 
by Weingartner & Draine (2003; hereafter WD03). They 
studied the alignment by RATs produced by a monochro- 
matic radiation field and for one particular radiation direc- 
tion ■(/) = 70°, taking into account thermal fiuctuations and 
thermal flipping. They observed the appearance of new at- 
tractor points at low angular momentum, but the underlying 
physics of this effect remained unclear. 

WD03 also posed a question concerning what is hap- 
pening when the entire spectrum of the interstellar radiation 
fleld (the ISRF) is accounted for, and when radiation arrives 
from other directions. Their study does not consider the ef- 
fect of gaseous bombardment and effects of H2 formation. 
We address these and other issues in the present below. 



* Although being a step forward compared to the earlier naive 
predictions of polarization which were mostly detached from the 
grain alignment theory, CL05 and the follow up studies (e.g. 
Pelkonen et al. 2007; Bethell et al. 2007; Cho & Lazarian 2007) 
are also not exact, as they are based on the alignment efficien- 
cies inferred from idealized numerical studies, rather than on the 
exact RAT alignment theory. 



The structure of the paper is as follows. In §2 we present 
current theoretical understandings of the RAT alignment 
and formulate our theoretical predictions. We briefly de- 
scribe thermal wobbling and our method of averaging RATs 
in § 3. In §4 we derive the analytical expressions for RATs for 
the AMO averaged over thermal fluctuations, and study the 
influence of thermal fluctuations to the grain alignment by 
RATs. In §5 we study the alignment for irregular grains with 
RATs from DDSCAT. We provide an explanation to the ap- 
pearance of the low attractor point based on the AMO and 
discuss the stability of low and high J attractor points in 
§6. We study the physics of crossovers induced by RATs in 
§7. We consider fast alignment in the presence of thermal 
fluctuations in §8. In §9 and §10 show the influences of ran- 
dom collisions and H2 formation torques in the framework 
of RAT alignment for both the AMO and irregular grains. 
An extended discussion is presented in §11. Finally, we sum- 
marize our results in §12. 



2 THEORETICAL CONSIDERATIONS 

In Paper I we provided analytical calculations for RATs for 
a grain model that consists of a perfectly reflecting oblate 
spheroid with a "massless" mirror attached at its pole (see 
the upper panel in Fig. [T} . This helical grain demonstrates 
RATs that are similar to those of irregular grains obtained 
by DDSCAT (see Paper I). 

The basic properties of RATs were the subject of a de- 
tailed analysis in Paper I. We established that the projection 
of RATs onto 63 axis (i.e. Qe3', see the upper panel in Fig.[l} 
determines the precession of the grain axis about the radia- 
tion direction. This component is also present for non-helical 
(e.g. simple spheroid) grains, and neither case induces grain 
alignment. It is the other two RAT components, Q^i and 
Qe2, that are responsible for grain alignment. In Paper I 
we explain why the RAT alignment tends to occur with the 
long axis of the grain either perpendicular to the magnetic 
field or perpendicular to the radiation direction (the choice 
of which, for a given external magnetic field, is determined 
by the ratio of the rate of radiative precession about the 
radiation direction and Larmor precession). 

In Paper I we found an important parameter that affects 
the grain alignment, the ratio QeT'"/Q™'', where 0^1""' 
and sae amplitudes of Qei and Qe2, respectively. Dif- 

ferent grain shapes illuminated by different radiation fields 
can have different ratios QS^'^/QS"'^, and thus the resulting 
alignment is different. 

We carried out the study of the RAT alignment in Pa- 
per I assuming that the grain always rotates about its prin- 
cipal inertia axis ai, corresponding to the maximal moment 
of inertia (hereafter called the maximal inertia axis). How- 
ever, this assumption is only valid for fast rotating grains. 
Such grains, are subjected to fast internal relaxation (Pur- 
cell 1979; Lazarian & Efroimsky 1999; Lazarian & Draine 
1999b) that couples the angular momentum J and ai . As the 
grain slows down to thermal angular velocity, the coupling 
becomes weaker. As a result, ai wobbles about J (Lazarian 
1994, Lazarian & Roberge 1997). Similar to DW97, in Paper 
I we disregarded this effect. 

Some implications of the thermal wobbling are self- 
evident. For instance, in Paper I we found that for a narrow 
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range of angles, when the radiation beam is close to be- 
ing perpendicular to magnetic field, the alignment becomes 
"wrong" (i.e. it occurs with the maximal inertia axis perpen- 
dicular to the mag netic field) B Such a "wrong" alignment is 
in contrast to what is widely observed in the ISM. However, 
we found that the "wrong" alignment in a diffuse medium 
corresponds to low angular momentum J. Therefore, we pre- 
dicted in Paper I that thermal wobbling of grains would 
destroy the "wrong" alignment. 

Qualitatively, in Paper I, we show that, in most cases, 
a high fraction of grains are aligned at a zero- J attractor 
point, assuming that ai is always parallel to J. In reality, 
thermal fluctuations induce the wobbling of ai with respect 
to J, resulting in the modification of RATs. This, in turn, 
changes the alignment of J with respect to B. Hence, we 
expect grains to be trapped at attractor points with ther- 
mal angular momentum J (hereafter called low-J attractor 
points) . 

One direct consequence of being trapped at the low-J 
attractor point is the decrease of the degree of the align- 
ment of grain axes with respect to the magnetic field. The 
degree of alignment is defined by the Rayleigh reduction 
factor (Greenberg 1968) 

i? = 1.5{cos^/3) - 0.5, (1) 

where (5 is the angle between the maximal inertia axis ai 
and the magnetic field B, (...) denotes the averaging over an 
ensemble of grains. Because the process of internal alignment 
between ai and J occurs much faster than the alignment of 
J with respect to B, _R could be approximately separated 
(see Lazarian 1994) as 

R^{R(e)){R{0). (2) 
Here {R{e)) and {R{C)) are given by 

(i?(6l)) = 1.5(cos^6i) -0.5, (3) 
{R{i)) = 1.5(cos^C) - 0.5, (4) 

where is the angle between ai and J, and ^ is the angle of 
J and B. 

When thermal fluctuations are absent, we have a per- 
fect alignment of ai with J, i.e. {R{6)) = 1. Therefore the 
Rayleigh reduction factor depends only on the degree of 
alignment of J with B. In the presence of thermal fluctu- 
ations, and assuming that the relaxation process obeys the 
normal distribution, (cos^S) is given by 

(cos^e) ^ Z cos'^e sin 6le-^'*''"='^''d6l, (5) 
Jo 

where Td is the dust temperature, Z = 
sin de~^^^^^'^'^'^d6 is the normalization factor, and 
E{e) = J^(l + {h- l)sm^e) /2hkTd is the kinetic energy 
of the spheroidal grain. Clearly, as the value of angular 
momentum J becomes comparable with the thermal value, 
Jth = \/1l\kTd where l\ is the maximal inertia moment 
of the grain, the aforementioned thermal wobbling should 
decrease the alignment of ai with B. 

On the basis of our findings in Paper I, we can qualita- 
tively address some questions posed in WD03. What will be 

^ "Wrong" alignment without specifying the conditions for it was 
also reported in DW97. 



the effect of random collisions on grain alignment? How do 
suprathermal torques, for instance, torques arising from H2 
formation (Purcell 1979) influence the alignment by RATs? 

We expect that grains trapped at low J attractor points, 
can be signiflcantly affected by random collisions of gas 
atoms. When the phase trajectory map of grains has attrac- 
tor points at high- J, we expect that random collisions can 
depopulate grains aligned at the lower attractor point, and 
stochastically direct grains to the high- J attractor point. As 
the grains with high J are immune to randomization and in- 
ternal alignment for high J is close to perfect (Purcell 1979), 
counter-intuitively, random gaseous bombardment can in- 
crease the degree of alignment. However, the removal of 
grains from the low-J attractor points is expected to make 
grain alignment time dependent, althoughh the correspond- 
ing time may be long compared with other processes in the 
system. With regards to H2 un-compensated torques (Pur- 
cell 1979; Lazarian 1995), because they are fixed in the grain 
body coordinate system, their effect depends on the fiipping 
rate of the grains. 

Obviously, at the high attractor point corresponding to 
J ^ Jth, the flipping rate is very low, therefore, H-2 torques 
act together with RATs and change the angular momen- 
tum depending upon the resurfacing process. In contrast, 
the grains flip fast in the process of heading to low attractor 
points. Assuming that the correlation time-scale is shorter 
than the alignment time, then grains are rapidly driven to 
low angular momentum at which the grain flips very fast; 
thus H2 torques will be averaged out to zero (see equation 
I62p . and could result in an additional randomization at low 
attractor point. 

Below we test our expectations with numerics. 



3 RATS AND EFFECT OF THERMAL 
WOBBLING 

In this section we briefly discuss the general definitions 
of RAT components, then present methods to average the 
torques over free motion and thermal wobbling. We also an- 
alyze the role of the third component Qes of RATs. 

3.1 RATs: spin-up, alignment and precession 

Similar to Paper I, in order to easily compare our results 
to those in earlier works, wherever possible, we preserve the 
notations adopted in DW97. 

As in Paper I, we consider only the anisotropic com- 
ponent of the radiation field, the radiative torque is then 
defined by 

r..,= "'""'^-^'^\ Qr(e,A3>), (6) 

where 7 is the degree of radiation anisotropy, a^jj is the 
effective size of the grain (see DW96; Paper I), and \,Urad 
are the mean wave length and mean energy density of the 
radiation field. The RAT efficiency vector Qr is represented 
in the lab system 61,62, and 63 

Qr(e, /3, $) = gei(e, $)ei ge2(e, $)e2 + Qeiie, P, $)e3, 

(7) 

where O, /3 and $ are angles describing the orientation of 
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Figure 1. Upper panel: Geometry for AMO consisting of a per- 
fectly reflecting spheroid and a weightless mirror. Lower panel: 
The orientation of a grain, described by three principal axes 
0.1,0.2,0.3, in the laboratory coordinate system (scattering refer- 
ence system) 61,62,63 is defined by three angles B,/9, The 
direction of incident photon beam k is along ei . 



ai in the lab system (see the lower panel in Fig. [1} . It was 
shown in Paper I that components Qei(0,/3, $ = 0) and 
Qe2{Q,P, "1> = 0) have universal properties and play a ma- 
jor role in the process of grain alignment, whereas the com- 
ponent Qe3(6,/3, "I> — 0) does not affect either spin-up or 
alignment, provided that ai is coupled with J. 

Here, we only deal with the alignment of angular mo- 
mentum with respect to the radiation direction or magnetic 
field (i.e. external alignment), therefore it is convenient to 
consider RATs in the spherical coordinate system J, ^, (j) (see 
Fig.H). 

In this coordinate system, the RAT can be written as 



lUraddef 



(8) 



where F, which is the torque component parallel to ^, acts 
to change the orientation of J in respect to B; H, the com- 
ponent parallel to J, is to spin grains up and G induces the 
precession of J about the magnetic field or radiation. These 



Figure 2. Alignment coordinate system in which tp is the angle 
between the magnetic field B and the radiation direction k, ^ is 
the angle between the angular momentum vector J and B, c/i is 
the Larmor precession angle. 



are given by 

F{iIj, 4>, S^) = (5ei(^, tp, (/))(— sin -(/jcos ^cos (p — sin ^cos ip) 
+ Qe2{^,ip,(p){cos -i/icos ^cos (p — sm ^sin ip) 
+ Qesi^, Ip, 0)cos ^sin (j), (9) 

G{Tp,(p,^) = Qei (Ci 'P)sm Tpsin (p — Qe2 {^,tp, (p)cos V'sin (p 

+ Qe3(C,V',0)cOS (10) 

H{ip, 4>, = <3ei(^, tp, 0)(cos ^cos ^ — sin tpsin ^cos (p) 
+ Qe2{i, tp, cp){smtpcos 5 -I- cos i/)sin ^cos (p) 
+ Qe3i^,tp,(p)sm ^sin (p, (11) 

where Qei{^,tp,(p),Qe2i^,tp,(p),Qe3{^,tp,(p), as functions of 
5, tp and (p, are components of the RAT efficiency vector in 
the lab coordinate system (see DW97; Paper I). To obtain 
Qei{(.,tp,<p),Qe2{^,tp,<p) and Qe3i^,tP,(P) from Qr(e,/3, $), 
which is calculated using the AMO or provided by DDSCAT 
for irregular grains, we need to use the relations between 
^, tp, (p and 9, /?, $ (see WD03 and Appendix A). 

3.2 Free-torque motion 

In the absence of external torques, a grain rotates around 
its principal axes. This motion is called free-torque motion. 
For a symmetric grain (e.g. the spheroidal body of AMO 
or brick; Paper I, Spitzer & McGlynn 1979), the free-torque 
motion consists of the nutation of angular velocity lu around 
J at a constant angle 7 and the rotation around the maximal 
inertia axis ai (see Fig. |3} with a period P-r- 

In general, an irregular grain can be characterized by 
an ellipsoid with moments of inertia 7i , I2 and I3 around 
three principal axes ai,a2,a3, respectively. For a freely ro- 
tating grain, its angular momentum is conserved (i.e., J is 
fixed in space) , while the angular velocity u nutates around 
J. We can call the wobbling associated with the irregular- 
ity in the grain shape irregular wobbling, to avoid confusion 
with thermal wobbling (also thermal fluctuations) induced 
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Figure 3. In body frame, the motion of a briclc consists of the 
nutation of f2 about ai and the precession of J around ai . 

by the Barnett effect (Purcell 1979) and nuclear relaxation 
(Lazarian & Draine 1999b). 

Obviously, the time-scales of the free-torque motion 
(e.g. rotation period P^) are much shorter than other time- 
scales (e.g., internal relaxations and the gas damping time; 
see Lazarian 2007). As a result, it is always feasible to aver- 
age RATs over the free-torque motion. A description of the 
free-torque motion for an asymmetric top in terms of Eu- 
ler angles 7, a and C, (see Fig. IA1|I can be found in classical 
textbooks (e.g. Landau & Lifshitz (1976; see also WD03). 
For an AMO with a spheroidal body, as in Fig. [1] the Eu- 
ler angle 7 is constant, while a and C, are functions of time 
(see Spitzer & McGlynn 1979). For irregular grains, we can 
get the Euler angles by numerically solving the equations of 
motion (see Appendix [C]). 

3.3 Averaging over free torque motion and role of 

As mentioned earlier, the rotation period Pt about the max- 
imal inertia axis and the precession time are both much 
shorter than the gas damping time. Therefore, we can aver- 
age RATs over these processes. 

In Figure |4] we show the phase trajectory of the tip of 
ai about J for a spheroid and a triaxial ellipsoid for two 
initial angles 70 = tt/IO and 70 = 7r/4. It is shown that in 
the case of spheroid, the trajectory is a circle corresponding 
to the precession of ai about J at a constant angle 7, but 
the evolution of the tip of vector ai produces a torus shape 
for the irregular grain. As a result, the average over free- 
torque motion for spheroid is concerned only the average 
over a circle, i.e. over the precession, whereas it is required 
to average RATs over the torus area (see Fig. |4ll . 

The averaging algorithm over free-torque motion for tri- 
axial elhpsoid has been realized in §3 in WD03. However, 
they averaged RATs over a single period of rotation P-r. 



irregular grain vs. spheroid 
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Figure 4. The evolution of the tip of the maximal inertia axis ai 
around the angular momentum for two initial angles 70 = n/10 
{upper panel) and 70 = n/A {lower panel) of ai and J. Here f 
is the nutation angle of ai about J and 7 is the azimuthal angle 
between ai and J. Filled dots and torus show the evolution of ai 
around J for the spheroid and the irregular grain respectively. 

We feel that a different averaging over a longer time is 
appropriate. This is well motivated due to a many orders of 
magnitude difference of the rotation time and the time-scale 
of internal relaxation (see Table [l|. 

We compare the results obtained by averaging RATs 
from an AMO (see equations (|23|l - (|25l )) with the body be- 
ing a triaxial ellipsoid over thermal ffuctuations for time 
step iV = 10^ and = 10^ with resuffs for Q^a = and 
A'^ — 10^ in Fig. \5\ It can be seen that the torque compo- 
nents {F)^, {H)<t> for the first and the last cases are almost 
the same, but their averaged torques have a small differ- 
ence with those obtained in the second case. Therefore, we 
can expect that the contribution of Qes to the spinning and 
ahgning torques is negligible when the averaging of RATs is 
performed with a sufficient accuracy (i.e. with sufficient high 
time step). In our study for the AMO, Qes is disregarded 
for both the alignment and spin-up process. 

3.4 Thermal fluctuations and thermal flipping 

Thermal fluctuations (i.e. thermal wobbling) and flipping 
arise from the coupling of rotational and vibrational degrees 
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AMO, ip = 70^ 
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Figure 5. (F)ij, and {H)^ are shown for different time-steps A'^ of 
averaging and for the case Qe3 = corresponding to two values 
of angular momentum. It can be seen that the torques obtained 
with N = 10^ steps do not depend on Qes. 



of freedom induced by internal relaxations. The effect was 
first discussed in Lazarian (1994). At that time, the strongest 
internal relaxation was believed to be associated with the 
Barnett effect (Purcell 1979). 

The Barnett relaxation can be easily understood. In- 
deed, a freely rotating paramagnetic grain acquires a mag- 
netic moment that is parallel to angular velocity fl due to 
the Barnett effect (first discussed in this context in Dolgi- 
nov & Mytrophanov 1976). The Barnett effect is the phe- 
nomenon of transferring the macroscopic angular momen- 
tum of a rotating body to electrons through flipping elec- 
tronic spins. To some degree, the magnetization of the ro- 
tating body is analogous with the body at rest, which is mag- 
netized by a rotating external magnetic field. The Barnett 
equivalent magnetic field is Hse = 5^/7 where 7 = e/2mc 
is the magnetogyric ratio of electrons (Purcell 1979). 

Because fl may not coincide with the maximal inertia 
axis ai, it precesses continuously in the grain coordinate 
system (see Fig. |3]). Therefore, the Barnett equivalent field 
can be decomposed into the constant component parallel 
to the precession axis and the rotating component which 



is perpendicular to 110. Apparently, the rotating component 
will induce a dissipation of rotational energy. As a result, we 
have an alignment of angular velocity and maximal inertia 
axis with angular momentum. 

Lazarian & Draine (1999a) found a new effect related 
to nuclear spins, which they termed nuclear relaxation. Sim- 
ilar to the rearranging of electronic spins, nuclear spins also 
become oriented by angular momentum transferred from 
the body. Although the nuclear magnetization arising from 
grain rotation is mostly negligible compared to that aris- 
ing from electrons (see Purcell 1979), the nuclear relaxation 
was shown to be much more efficient than Barnett relaxation 
for 10~^ < a^ff < 10^^ cm grains. This can be easily un- 
derstood. Spin flipping is a mechanical effect that depends 
on the angular momentum of the species rather than on 
their magneto-magnetogyric ratio 7. Thus, a rotating grain 
will have the same portion of nuclear and electron spins 
flipped. The magnetic field that induces such a fiip is differ- 
ent for electrons and nuclei (i.e. it is inversely proportional 
to 7). The dissipation is proportional to the "equivalent" 
field squared (i.e. to B^^ ~ -ij). It does not depend on 
the value of the nuclear magnetic moment, as the lag be- 
tween magnetization and fl increases with the decrease of 
the magnetic moments. In other words, the coupling between 
nuclear spins is less efficient than the coupling between elec- 
tron spins; hence, there is a substantial lag in the nuclear 
spin alignment when f2 precesses around ai. All in all, al- 
though the magnetic moment arising from nuclear spins of 
a rotating body is negligible, the corresponding relaxation 
(i.e. nuclear relaxation) is approximately 10^ times higher 
than the Barnett one. 

Thermal fluctuations within the grain body (see Purcell 
1979) coupled via internal relaxation with the macroscopic 
rotation of the grain can affect the internal alignment, and 
result in random deviations of the major axis ai in respect 
to J. Following the fluctuation-dissipation theorem (Landau 
& Lifshitz 1976), the thermal equilibrium distribution of ai 
deviations can be established. The average of a torque A 
over thermal fluctuations for a spheroid with Ii > I2 ~ I3 
is defined by (Lazarian & Roberge 1997) 



/;/'d7A(7, J) sin -ye^-E/kTa] 



{A} = 



Ji^^^ dj sin 7exp[— i5/fcTd] 



(12) 



where h — j^, 7 is the deviation angle of ai and J, and the 
kinetic energy E of the grain is 



i?(7) 



2/1 



[1-f (/i-l)sin'7]. 



(13) 



For irregular grains with Ii > I2 > I3, the corresponding 
average is given in Appendix [Cl 

Using the RAT expressions obtained for the AMO (see 
Paper I), one can explicitly integrate equation (|12p to ob- 
tain RATs induced by thermal fiuctuations. However, for 
irregular grains, we need to numerically average RATs ac- 
cording to equation (Dl). The resulting averaged RAT 
components are plugged into equations (|9llip to obtain 



^ A more accurate description of the process that account for 
the finite spin-lattice relaxation is provided in Lazarian & Draine 
1999b 
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(B = 5 ij,G) 
(B = 50 /iG) 




(J/J.J' 



Figure 6. Larmor precession, internal relaxation and thermal 
flipping time-scales for B = 5fj,G and B = 50fiG for a grain size 
of file// = 0.2/im. When magnetic field increases, the Larmor 
precession time-scale decreases and becomes comparable to the 
Barnet relaxation time-scale for small J. 



{F{^,<l>,%b,J)),{H{^,<P,^,J)},{G{^,<l>,^P,J)), which are re- 
quired to solve the equations of motion. 

The internal relaxation can induce the axis ai to flip 
over with respect to J. This phenomenon is called ther- 
mal flipping (Lazarian & Draine f999a). The probability 
of thermal flipping has been obtained in Lazarian & Draine 
(f999a), 



tr/ = t-iexp[-0.5(^ - f)] 



(14) 



where Jth = \/2/i kTd is the thermal angular momentum 
corresponding to the dust temperature, tsn = (^ijar + 
^riuci)~^ is the internal relaxation time for Barnett and nu- 
clear relaxation processes. Because t„uci is usually much 
shorter than tsar for grains from 10"'^ to 10"" cm, the nu- 
clear relaxation dominates the process of thermal flipping for 
typical aligned interstellar grains (see Kim & Martin 1995). 

In our analysis, we assume that thermal fluctuation and 
flipping time-scales are shorter than the Larmor precession 
time. Fortunately, for the magnetic field of the ISM and for 
astronomical silicate material, the Larmor precession time- 
scale is much longer than that of thermal fluctuations and 
thermal flipping. Therefore, averaging over the later motion 
is appropriate. 

For circumstances in which magnetic field is stronger, 
the Larmor precession time-scale may be comparable with 
the thermal fiipping time-scale. For instance. Fig. [6] shows 
that for B = bOfiG and for ordinary paramagnetic grains, 
tLar ^ ttf for J < lOJth- FoT this case, in order to treat 
properly the grain dynamics, it is necessary to follow both 
the thermal fluctuation and Larmor precession. 



3.5 Equations of motion 

The motion of a grain subjected to a net torque is completely 
determined by three variables: the angle ^ between the an- 
gular momentum vector J and the magnetic field direction 
B, the precession angle of J around B and the value of the 



angular momentum J (see the lower panel of Fig. 1). The 
equations of motion, if we disregard paramagnetic dissipation 
for these variables are 



^"rad^eff'^ 
2Jsin ^ 



dt 



d^ ^ 7"rad"eff-^ 
2J 



{F(C,0,^,J)), 



dt 

= ^7%ad«eff^(-'^ "^)) 



(15) 

(16) 
(17) 



where {F(C, V, J)), {G(e, 0, ^Z-, J)) , and {H{^A,^,J)) are 
defined by equations (f9)l- (|ll|l . which are already averaged 
over thermal fluctuations (see WD03) , tgas is the gas damp- 
ing time-scale (see Table[l| and f2_B is the Larmor precession 
frequency of the angular momentum about the magnetic 
field. 

For the ISM, the Larmor precession time is always 
shorter than the gas damping; thus, we can average equa- 
tions (|16|) - (|17|) over a precession period. As a result, equa- 
tions p5p - p7p can be reduced to a couple equation for ^ and 
J only in which the spinning and aligning torque (F) and 
(H) are replaced by {F)if, and {H)^, which denote the quan- 
tities obtained from averaging corresponding RATs over the 
Larmor precession ang le 4, from to 27rl!| 

A stationary point in a phase trajectory map is deter- 
mined by ^s, Js, which are solutions of the equations of 
motion (see DW97) 



d^ 
dt 
dJ_ 
'dt 



= 0, 
= 0. 



The above stationary point is an attractor point if 



1 d{F}^ 



(H)^ d^ 



<0, 



(18) 
(19) 

(20) 



and is a repellor point otherwise (see DW97). 

From equations p6p and <\17^ . we get ^s, Js that satisfy 



{F{^s,i,,J))^ = 0, 



(21) 



Js = tgasT^JUj.g^^algX{H{£,s,'ip,J))4.. (22) 



4 RAT ALIGNMENT FOR AMO 

We first consider the role of thermal fluctuations in grain 
alignment based on an AMO consisting of a reflecting 
spheroid {I2 ~ I3, hereafter spheroidal AMO) and a mir- 
ror as in Paper I. Then, in order to see the correspondence 
of the AMO with irregular grains in terms of dynamics, we 
replace the spheroid by an ellipsoid with the principal mo- 
ments of inertia Ii > I2 > I3 (hereafter ellipsoidal AMO). 



^ For the sake of simplicity, hereafter, we denote 
(F) = (F(C,</.,V,J)), {H) = {H{^,,i>,^,J)), and 
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4.1 RATs: general expressions 



4-2.1 Analytical averaging RATs for one component Qei 



In terms of RATs, an AMO is formally only applicable for 
A <C Uef f as only in this case, is geometric optics approach, 
used to derive the analytical formulae in Paper I, appropri- 
ate. However, in Paper I we proved that the functional forms 
of RATs for the AMO and irregular grains for A ^ Oef / are 
similar. By choosing the appropriate ratio Q^^^/QS"^ it is 
possible to see that the dynamics of the AMO is similar to 
that of irregular grains, if their torque ratio Q^^'^/QS"'^ is 
the same. Therefore, the AMO can act as a proxy for ac- 
tual grains, and the advantage of this is that it allows an 
analytical insight into grain alignment. 

For an AMO in which the mirror is tilted by an angle 
Q with the axis 3.2, the RAT components are given by (see 
Paper I) 



Oei(e,/?,o) 



-—^{\ni cos O — n2 sin O cos P\[nin2 cos 
A 

2 2 
- cos 13 sin 20 ^ cos (3 sin 20 



2 2 
nin2 sin ^0 cos ^/9]), 
4/1 



(23) 



Qe2(6,/5,0) = — (|ni cos — 712 sin O cos f3\[ni cos f3 cos 
A 



"1^2 2r, ■ r,r, 

cos p sm 20 

2 ^ 

-I- cos /3 sin ^0]), 



nin2 



sin 20 



(24) 



4li 

Qe'iiQ, P, 0) = —r-(\ni cos — 712 sin cos f3\ni sin /3 
A 



X [ni cos — 772 cos (3 sin 0]) 

l)K{e) sin 20, 



, fe , 2 2ea , 2 



A 



(25) 



where is the angle between the axis of major inertia ai and 
the radiation direction k (see Fig. [BoTTjer); h is the distance 
from the square mirror of side h to the center mass, A is 
the wavelength, 77i = — sin a, 772 = cos a are components 
of the normal vector of the mirror in the grain coordinate 
system (see Fig. [T][f|, a, b are minor and major semi-axes of 
the spheroid, s — a/b < 1 and e is the eccentricity of the 
spheroid; K{Q) is the fitting function (see also Paper I). As 
in Paper I, we treat the AMO with a — 45° as our default 
model. 

RATs at a precession angle <1> (see Fig.[T|ow)er) are given 
by (see DW97) 

Qei(0,/?,-l>) = Qei(0,/3,O), (26) 
Qe2(0,/3,$) =Q.2(0,/3,O)cos$-Q,3(0,/3,O)sin-l', (27) 
Qe3(0,/3,$) = Q.2(0,/3,O)sin$-f Q,3(0,/3,O)cos<E>. (28) 



4.2 Alignment with respect to k 

First we consider the alignment in respect to the direction 
of radiation k, which also correspond to the situation when 
the direction of light k coincides with that of magnetic field 
B (i.e. tp = 0°). 



As discussed earlier in §3.4, Qe3 does not affect RAT align- 
ment, apart from inducing the precession of angular momen- 
tum J about k. Thus, the alignment problem only involves 
Qei and Qe2- To clarify the role of the torque components, 
we first consider the alignment driven by the component Qei 
in the presence of thermal fluctuations. 

For the default AMO (i.e. a = 45°), the contribution 
arising from the change in the cross-section is insignificant 
(see Fig. 6 in Paper I). Thus, we can ignore the factor A± = 
|ni cos — 772 sin cos /3j present in equations (|23|) - (|25p . 
As a result, with an accuracy of 5%, equation (|23|l can be 
rewritten as 



Qci ~ Q\ 



•iaxr(3cOS^0 - 1) 



COS COS (3 -I- sin cos 2/3] , 
(29) 



where QS"" = ^nin2. 

When there is incomplete alignment of ai and J, we 

have 



cos = cos ^ cos 7 — sin ^ sin 7 sin 77, 



(30) 



where 77 is the precession angle of ai around J, 7 is the angle 
of ai and J, and ^ is the angle between J and B. In addition, 
13 is & complicated function of Euler angles. 

Substituting equation pop into (|29p and averaging the 
resulting expression over the precession angle r], the second 
and third term involved (3 are averaged to zero. Hence, we 
obtain 



^ ^ r^max,ni 2 2, , siu 7sin ^ 

Qei ~ Qei [3(C0S 7COS ^ H 



)-!]■ (31) 



Using equations '\V2\ and (|13p for the expression (|3ip we 
obtain 



{Qei) ~ 



j;d7(cos%cos^C+ ''"'f"'^ 



-1) 



-^[i+{/i-i)sin 7] 
J(j (17 sm 76 th 



X sm 76 



-[i+{h-i)sm 7] 



(32) 



where Jth ~ \/2/i kTd- The aligning and spin-up torque are 
respectively given by (see Eqs. 111)113^ 



° Note that RATs in equations I I23I25I I have opposite signs com- {F)4> — ^(Qei) sin ^ (33) 

pared with those in Paper I because in Paper I we have defined 
ni = sin a (i.e. we incorporated the minus sign of RATs into ni) 



{H)4, = {Qei) cos e (34) 
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Hence, 



{F)^,^- smiQ^i — 



rze -^f'l (1 + 3 cos 20fe th yJ^-V^Th^l 

lax V ^ L V ^ 



+ sm ^y^i 



(3 + 2(-l + h)j^)]^Erfi{V-l + ftj^) ^ 



(35) 

,2 ,2 A 
(-l+h)^ hAy- , , ^ 

r= e -'th (1 + 3 cos 20 [e V-l + hi^ v 

[n)^ ^ COS ^y^i — 



J§-(-l + h)e~ ^Erfi{^-l + hjL.) 



e^(3 + 2(-l + h)-l^)]y^Erfi{V~l + 
cosCQT" 



J 2 



-(-l + ft)e ■'*.i5r,fi(^/-l + /i^) 

(36) 



Equation (|35|l reveals that {F)^ = for ^ = 0, vr and 
cos 2^ = —1/3, regardless of angular momentum J. In ad- 
dition, zero points cos ^ = and cos 2^ = —1/3 of (i?),^ 
do not depend on J either. This indicates that thermal fluc- 
tuations within the spheroidal grains do not alter the value 
of angular momentum (i.e. J = 0) at low-J attractor points 
produced by RATs when J — ^ Jth- Therefore, we expect 
that the resulting alignment is not significantly affected by 
thermal fluctuations. 

The above features of RATs can be seen in Fig. [T] 
which shows RAT components {F)^ (solid line) and (i?),^ 
(dashed line) as functions of ^ at a value of angular momen- 
tum J = Q.lJth- It follows that there will be four station- 
ary points in the phase map, corresponding to {F)^ = 
at 5 = 0, 7r/3, 27r/3, TT. In addition, the stationary point 

f = is a high- J attractor point as /if ^t. < and 

(77)^(5 — 0) > (see the upper panel in Fig. [7]). The lower 
panel of Fig. [7] shows (F)^ and (H)^ as functions of J/ Jth 
for ^ — n/3. There, it can be seen that for J ^ Jth, {F)^ 
and {H)4, are saturated as a result of the perfect coupling 
of ai with J, which makes RATs independent on the angu- 
lar momentum. As J j Jth decreases, (F)^ and {H),^ decrease 
steeply (see Fig. [Tjl . 

^.2.2 Averaged RATs for both components Qei,Qe2 

The analytical averaging over thermal fluctuations for Qe2 is 
more complicated because of its dependence upon which 
is a function of Euler angles a, 7 and f , <j), V" (see equation 
[57] and Appendix A). Therefore, we use numerical averaging 
for RATs, rather than deriving analytical expressions for 
them. The resulting torques {F),p, (H)^ are used to solve 
the equations of motion p6|l - (|17p . 

Fig. [8] shows spin-up and alignment torques for J = 
10""^ Jth, 0.9 Jth and lOJth corresponding to cases in which 
thermal fluctuations are dominant, important and negligi- 
ble, respectively. It can be seen that for J 3> Jth, {F)^, — 
for cos £, = =f1. It turns out that for J 3> Jth, the AMO 
creates only two stationary points corresponding to per- 
fect alignment of J with B. When J — 10~^Jth, there are 
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Figure 7. Spheroidal AMO: aligning (F)^, and spinning {H)^ 
torques averaged over thermal fluctuations for the case ip = 0° 
and = 0. Upper panel shows {F)^ and {H)^ as functions 

of t; for J = O.lJth- Lower panel shows that (F)^ and (H)^ for 
^ = 7r/3, decrease rapidly with J/ Jth decreasing and get saturated 
as J/ Jth > 1- 

new stationary points at cos ^ = ±0.9 corresponding to 
(F)^ = 0. This means that thermal fluctuations can produce 
new stationary points. At the same time, the magnitude of 
spinning torque {H}cf, decreases as J decreases (i.e. when 
thermal fluctuations increase; see Fig. [Hfouier). As a result, 
the alignment of grains is expected to be similar to the case 
without thermal fluctuations. 

4-2.3 Trajectory maps 

Let us consider the dynamics of the spheroidal AMO driven 
by RATs averaged over thermal fluctuations, and test the 
predictions using the analytical results above. 

We solve the equations of motion for grains having the 
same initial angular momentum and uniform orientation dis- 
tribution with respect to the interstellar magnetic field. Pa- 
rameters necessary for calculations for the ISM conditions 
are given in Table [1] Phase trajectory maps are plotted in 
coordinates ( cos (^,J / I\lot) where ^ is the angle between J 
and B, and ujt ~ ^fWIfTx. The upper and lower parts of 
which correspond to ai initially parallel and anti-parallel to 
J, respectively. For a grain size a^// = 0.2iJ,m and shape 1 
(see Fig. llSp , the ISRF can produce high stationary points 
with Jhigh ~ 200/1 tJT for ?/> = 0°. To represent high and 
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Table 1. Physical parameters for diffuse ISM 
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Figure 8. Spheroidal AMO: Aligning [upper panel) and spinning 
(lower panel) torques averaged over thermal fluctuations for dif- 
ferent J/Jti^. Upper panel shows that as J = 10~^Jf^, one new 
zero point of (i^)0 appears at cos ^ = —0.9. Both (F)^ and (H)^ 
exhibit rapid decrease with J/Jth decreasing. 



low attractor points together, we consider the ISRF with 
Urad ~ uisrf/W where uisrf = 8.64 x 10~^^erg/cTn^ is 
the energy density of the ISRF. Thus, the phase trajectory 
maps are shown with Jmax = 20IiU)t in the present paper. 

Further in the paper, a stationary point on the phase 
trajectory maps is marked by a circle is an attractor point, 
which is the point to which adjacent trajectories tend to 
converge, and a stationary point marked by a cross denotes 
a repellor point (i.e., the trajectories are repulsed while ap- 
proaching it; see also Paper I)Q In general, a phase trajec- 
tory map may have low-J and high-J attractor points or 
only low-J attractor points (see more in Paper I). To see 
clearly the modification induced by thermal fluctuations on 
grain dynamics, we frequently show side by side the tra- 
jectory maps for the case without thermal fluctuations as 
in Paper I, and with thermal fluctuations. Phase trajectory 
maps for the spheroidal AMO are shown in Figs. |9] 1101 and 
111! in which the upper and lower panel correspond to the 
cases without and with thermal fluctuations. 



^ In some trajectory maps, we do not label all existing repellor 
points. 



Definitions 



Values 



Gas density 
Gas temperature 
Gas damping time 
Dust temperature 
Anisotropy degree 
Mean wavelength 



n = 30 cm ^ 



tgas =4.6 X 10l2(_-£-^)a_5 



Td = 20 K 
7 = 0.1 
A = 1.2 /im 

Mean density of thelSRF uisRp = 8.64 X 10^^^ erg cm"^ 
Effective grain size a,^f{ = 0.2/im 



Here Tgas = Tgas/WO K, n = n/30 g cm^^, and a_5 = 
a^f f /10~^ cm. p = p/3 g cm~^ where p = 3 g cm~^ is the 
mass density of the grain. 



AM0,'^=0°,Qg2='-'> thermal fluctuations 




AMO, Q,2=0, ^ = 0°. thermal fluctuations 




Figure 9. Spheroidal AMO: Phase trajectory maps for the align- 
ment by Qei only. The upper panel shows the phase map in 
the absence of thermal fluctuations and the lower panel shows 
the phase map when thermal fluctuations are accounted for. The 
panels show that the position of the zero-J attractor point at 
cos g = —0.6 is not changed even in the presence of thermal 
fluctuations. 
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AMO, -^ = 0°, thermal fluctuations 




-1.0 -0.5 0.0 0.5 1.0 



Figure 10. Spheroidal AMO: similar to Fig. [9] but grains are 
aligned by all torque components. Figs show the shift of the zero 
J attractor point from cos § = —1 {upper) to cos J = —0.9 
{lower), but its angular momentum remains equal zero even when 
thermal fluctuations are accounted for. 

When the grain ahgnment is only driven by Qei and 
tp = (f (Fig.H]), we see that each phase map has two high- J 
attractor points A and B. In the absence of thermal fluctu- 
ations, the torque {H)^ decelerates grains to the attractor 
point C with J = (upper panel). When thermal fluctua- 
tions are present, the attractor point C is not affected. This 
result is consistent with our analytical predictions in §4.2.1 

When the components Qei and Qe2 act together. Fig. 
1 101 shows that the angular momentum of low attractor point 
remain the same (i.e. J remains equal zero). However, its 
position is slightly shifted to cos ^ = ±0.9. 

4.3 Alignment with respect to B 

Below we consider grain dynamics when the magnetic held 
plays the role of alignment axis. As an example, the radiation 
direction t/j = 70° is adopted. 

Fig. [11] shows that thermal fluctuations do not increase 
the value of angular momentum at the attractor point C for 
the case t/j = 70°. For other angles xp, we also found that 
the zero-J attractor point C is unchanged in the presence of 
thermal fluctuations. 

It is easy to see that in the assumption of ai||J, the 



AMO, '^^ = 70°, no thermal fluctuations 




-1.0 -0.5 0.0 0.5 1.0 



AMO, "^=70°, thermal fluctuotions 




-1.0 -0.5 0.0 0.5 1.0 



Figure 11. Spheroidal AMO: similar to Fig.[TOl but for ^ = 70°. 
Figs show that the zero J attractor point C at cos ^ = —0.9 
(upper panel) panel) is unchanged in the presence of thermal 
wobbling (lower panel). 

irregular shape of the ellipsoid of inertia (i.e. /i 7^ /2 7^ 
Is), is not important for the grain dynamics. However, the 
irregularity in the grain shape becomes important in the 
case of a wobbling grain because the averaging of RATs 
over free-torque motion (see Fig. |4| depends on its ellipsoid 
of inertia. To address such effects, the spheroidal AMO can 
be modifled. For instance, if the mirror is not weightless, 
the free motion of the spheroidal AMO will be of a triaxial 
ellipsoid of inertia, rather than a spheroid. Further, in ij4.4l 
we replace the spheroid of the AMO by a triaxial ellipsoid. 



4.4 Alignment for an ellipsoidal AMO 

In this section, we study the alignment of the ellipsoidal 
AMO in which the spheroid body (see Fig. [T} is replaced by 
an ellipsoid with moments of inertia Ii : I2 ■ H = 1.745 : 
1.62 : 0.876, which is similar to shape 1 (see Fig. I15|l . 

We recall that stationary points, which determine the 
alignment of grains correspond to {F)^ = (see equa- 
tion [22]). Fig. 1121 shows the torque components for ijj — 70°. 
It can be seen that for J Jth (i.e., thermal fluctuations are 
negligible), {F}^ = for cos ^ = =Fl. This indicates that 
the ellipsoidal AMO produce two stationary points corre- 
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AMO, iy = 70°, right helicity 



AMO, -^^=70°, thermal fluctuations 
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Figure 12. Ellipsoidal AMO: aligning and spinning torques av- 
eraged over thermal fluctuations as a function of ^ for several 
angular momenta. Upper panel shows that as J = IJth^ one new 
zero point of (-F)^ appears at cos ^ = —0.8. While the change 
in sign of (H)^ in the vicinity of cos ^ = —1 for J = O.lJth is 
shown in lower panel. 



sponding to perfect alignment of J with B. As J ^ O.lJth, 
there appears a new stationary point at cos ^ = —0.8. It 
means that thermal fluctuations produce a new stationary 
point at cos ^ ~ —0.8. However, this new stationary point 



is a repellor point as 



COS es=-0.8 



> (see Fig. fT^ . 



{H}4,di I 

In addition, the lower panel shows the change in sign of 
the spinning torque from negative to positive as J ~ O.lJth 
for cos ^ ~ —1. Thus, the stationary point cos ^ — —1 is 



an attractor because 



< (see Fig. 



Therefore, thermal fluctuations within the ellipsoidal AMO 
modify the properties of RATs as J — > 0, and split a zero- 
J attractor point in absence of thermal wobbling into two 
low-J attractor points. The corresponding trajectory maps 
for this case are shown in Fig. [13] 

Now we rescale the amplitude of Qei , Qe2 so that 
QS^^'/QS""' = 0.78. For this AMO, Paper I shows that the 
grains are aligned with two high-J attractor points and a 
zero- J attractor point. However, in presence of thermal fluc- 
tuations and for the ellipsoidal AMO with I2 I3, the zero- J 
attractor point becomes the attractor points with the value 



Figure 13. Ellipsoidal AMO: the phase map for ip = 70° shows 
two low attractor points C and D' with J = 0.5/iu)t and two 
rcpcUor points A and B. 



AMO, -^^-70°, thermal fluctuations 




Figure 14. Ellipsoidal AMO: the phase map when the torque 
ratio QJJ5^"^/Q™^ is changed to the default value to 0.78 for 
^ = 70° shows two low attractor points with J = 0.5Jth ^'rid two 
high J attractor points A and B with J = 18 Jth ■ 



of angular momentum of the order of thermal value, as seen 
in Fig. [H 



5 RAT ALIGNMENT FOR IRREGULAR 
GRAINS 

5.1 Properties of averaged RATs 

We use the publicly available DDSCAT code (Draine & 
Flatau 2004) to calculate RATs for irregular grains (i.e. 
shape 1 and 3; see Fig. llSp . The optical constant function 
for astronomical silicate is adopted (Draine & Lee 1984). We 
computed RAT efficiency Qa(0,/3, "3?) in the lab coordinate 
system (see Fig. [T|, over 32 directions of O in the range 
[0, n] and 33 values of (3 between and 27r for $ = (an- 
gle produced by ai and the plane 61,63 where are three 
unit vectors of the lab coordinate system; see Fig. [T|. The 
RAT efficiency Qa(0,/3, '&) for an arbitrary angle "I> is eas- 
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DDSCAT, X=1.2/i,m, 0^,^ = 0. 2/i.m 
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Figure 15. Irregular grains of shape 1 and siiape 3 (see DW97) 
are taken to study RAT alignment by DDSCAT. 



ily calculated using equations (|26|l - (|28|l . In our paper, we 
calculate RATs for irregular grains of size a^ff = 0.2^m in- 
duced by monochromatic radiation field of A = 1.2/im and 
the spectrum of the ISRF (see Mathis et al. 1983). 

Fig. 1 161 show {F)^ and {H)^ (i.e., aligning and spinning 
torques) obtained by averaging the corresponding expres- 
sions (see equations [9l and 1 1 1 |l over thermal fluctuations for 
different angular momenta, and for the monochromatic ra- 
diation at angle il> = 70° toward B. 

From Fig. [16] we see that for J/ Jth 3> 1, = at 

cos^ = ±1, —0.65, corresponding to three stationary points. 
However, when J/ Jth decreases (i.e., thermal fluctuations 
become stronger), the intermediate stationary point shifts 
to the left, and disappears as J/ Jth —> 0. 

Fig. [17] is similar to Fig. 1161 but presents torques for 
the entire spectrum of the ISRF and -0 = 0°. It is also 
seen that the intermediate stationary point (i.e. point with 
{F}^ — 0) shifts to the left as J/ Jth decreases. It disappears 
when J /Jth = 0. Furthermore, following both Fig. 1161 and 
1171 (lower panels), we see that the spinning torque {H},/, de- 
creases with J decreasing. This is because for low angular 
momentum, thermal fiuctuations become stronger, so the 
axis of grain ai fluctuates with a wider amplitude around J. 



Figure 16. (i^)^ and (^^)<^ for the monochromatic radiation field 
at ■)/> = 70° as functions of the angle ^ between J and B for three 
values of J. 



As a results, RATs decrease (see analytical results for the 
spheroidal AMO in the upper panel of Fig. [7]l . 

Now let us consider the property of the intermediate 
zero point of the aligning torque (F)^ shown in the lower 
panel of Fig. 1171 We can check that this stationary point 
corresponding to J/ Jth ~ 2.5 is an attractor point. Accord- 
ing to equation (|20|) . the criteria in order for one stationary 



point to be an attractor point is 



< 1. With the 



stationary point at cos^s — —0.85, we have {H)if, > (see 
the upper panel), and < due to the decrease of 

(F)^ with ^ in the vicinity of (see the upper panel of Fig. 
1171) . Thus, TST-^ < 0, which satisfies the criteria for an 
attractor point. 



5.2 Phase Trajectory Maps for shape 1 

Let us consider the trajectory map for grains of size a^f / = 
0.2/im driven by RATs produced by a monochromatic ra- 
diation of A = 1.2^m, in the direction t/i = 70°, which is 
similar to the setting in WD03. As we mentioned above, an 
important difference in our treatment and that in WD03 is 
that we average RATs over 10^ rotation periods rather than 
over only one period as WD03 did. As we discussed in §3.3, 
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shape 1, X=^.2f^^r\, no thermal fluctuatior 
20 



DDSCAT, ISRF, a^ff=0.2/im 



0.010 




1 .0 



Figure 17. (F)^ and (H)^ for the ISRF at i/i = 0° as functions 
of the angle ^ between J and B for three values of J. 



one of the consequences of this is that the contribution of 
Qe3 component vanishes for ahgning and spinning torque 
components. 

The upper panel in Fig. [18] shows the trajectory map 
for the case without thermal fluctuations with two high- 
J attractor points (A, B) and one zero-J attractor point 
C. When thermal fluctuations are accounted for, the lower 
panel reveals the split of the attractor point C to C and 
D' as also seen in WD03. However, the value of angular 
momentum at C and D' is J = 0.7Jth, which is lower that 
the value obtained by WD03. This difference stems from the 
fact that the contribution of Qe3 to the spinning and aligning 
torques is completely canceled when we average RATs over 
a sufficiently long time. In fact, we tested that, if Qes is set 
to zero, the resulting map is the same with the lower panel 
of Fig. 1181 Yet, if we adopt the same averaging as in WD03, 
we also get the low attractor point with similar value J as 
in WD03. We also see a correspondence between our results 
here with those obtained with the AMO (see §4.4). 

The lower panel in Fig.[T8]also shows that about 20% of 
grains are aligned at two high-J attractor points A, B, and 
about 80% grains are driven by RATs to the low- J attractor 
points C and D. There, grains flip to the opposite flipping 
state (i.e. from upper to lower frame) and back. However, 
immediately after entering the opposite flipping state, grains 



+20 




+20 



+ 10 



shape 1, X=1.2 fJ-m, thermal fluctuation; 




+20 



Figure 18. Phase trajectory maps for A = 1.2/im and a = 0.2/jm, 
and ip = 70° corresponding to no thermal fluctuations {upper 
panel) and with thermal fluctuations {lower panel). For the case in 
which J is fixed, when grains are driven to the crossover J = 0, ai 
flips over. So grains change frequently their states corresponding 
to switching between the upper and lower frame in the map. 



are decelerated by RATs to that point again and they flip 
back. So grains flip back and forth between C and D. In fact, 
C and D are indistinguishable because grains flip very fast 
between them with the probability of finding grains on each 
point is ~ 0.5. In that sense, the points C and D are also 
crossover points 

Figs [19] and [20] show the phase maps of grain alignment 
for the full spectrum of the ISRF instead of monochromatic 
radiation, and for two directions of radiation — 0° and 90°. 
For the direction tp = 0° , the map in the lower panel consists 
of three attractor points A', B', C in which the point C is 
the old attractor point at very high angular momentum, and 
A' , B' are low- J attractor points originated from A and B at 
J = for the case of no thermal fiuctuations (see the upper 
panel). In other words, averaging over thermal fluctuations 



These crossover points are different from those described in 
DW97, in which the grains were spinning up right after crossover. 
This stems from inaccurate treatment of crossovers in DW97, see 
more detail in Paper I. 
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ISRF, DDSCAT, shape 1, B = 
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Figure 19. Trajectory maps for the entire spectrum of tfie ISRF 
and irregular grain of size a^f f = 0.2/jm when thermal fluctua- 
tions are absent (upper panel) and when thermal fluctuations and 
thermal flipping are included (lower panel), corresponding to the 
alignment with respect to k or i/) = 0° (equivalent to the grain 
alignment in the absence of magnetic field). 
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Figure 20. Similar to Fig. 1191 but for ■)/' = 90°. Upper panel: the 
map with two high J attractor point A, B, two zero-J attractor 
points C, E and a "wrong" attractor point D at cos 5 = 0. Lower 
panel: with thermal fluctuations and the "wrong" attractor point 
D disappears. 



of grains for the ISRF and = 0° also changes the zero-J 
attractor point to the new thermal J attractor point. 

In Paper I we pointed out that at tp close to 90° for 
both the AMO and irregular grains, the grain alignment 
tends to take place with ai _L B, which is in contrast to the 
classical Davis-Greenstein (1951) expectations. We termed 
this "wrong" alignment. An attractor point D at "wrong" 
alignment angle (i.e. ^ = 90°) shown in Fig, \20h vver for the 
direction ip = 90° together with four other attractor points 
A, B, C and E, in the absence of thermal wobbling (see also 
Paper I). However, in the presence of thermal wobbling, the 
"wrong" attractor point D is no longer existing. Instead, we 
obtain four new attractor points A', B', C and D', corre- 
sponding to J ~ Jth and cos ^ = ±1 (see the lower panel 
in Fig. I20p . It means that the "wrong" alignment is indeed 
eliminated by averaging induced by thermal wobbling (see 
also §2). 

5.3 Dynamics of shape 3 

For the sake of completeness, let us study the effect of ther- 
mal fluctuations on RAT alignment for another irregular 



grain (shape 3 in DW97). As an example, we consider a 
particular grain size a^ff = 0.2/im and one light direction 
ip = 70°. 

Fig. [21] shows the trajectory maps obtained for this 
grain corresponding to the case without (upper panel) and 
with thermal fluctuations (lower panel) . It is shown that for 
tp = 70°, this grain produces a phase map with two high- J 
attractor points A, B and two low- J attractor points C and 
D' in the presence of thermal fluctuations (see lower panel 
of Fig. I2ip . The lifting of the low- J attractor point C (see 
the upper panel of Fig. from J = to J = 2Jth (C, D') 
when thermal fluctuations are taken into account, is also 
seen for this shape. However, grains can rapidly flip back 
and forth between C and D'. 



6 LOW J ATTRACTOR POINTS 

Our study indicates that a high fraction of grains are aligned 
at low- J attractor points. Therefore, the origin and stability 
of this low-J attractor point are important for grain align- 
ment. We address these questions below. 
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Figure 21. Phase maps for the grain shape 3 and the ISRF 
with two high- J attractor points corresponding to the cases with- 
out thermal fluctuations (upper panel) and with thermal fluctua- 
tions {lower panel). Thermal fluctuations shift the zero-J attrac- 
tor point C to C, D' with J = 2Jth = 2Iiuit- 

6.1 Origin 

The results above show that thermal fluctuations produce a 
new low-J attractor point from the zero- J attractor point, 
for both the ellipsoidal AMO and irregular grains, but this 
effect has not been seen for the spheroidal AMO. Now let 
us explain why this occurs by using the AMOs, but the 
explanation can be applicable to irregular grains. 

First of all, let us study RATs averaged over the free- 
torque motion for the spheroidal and ellipsoidal AMOs. Fol- 
lowing equation (IDlfl (see Appendix D), the average over 
thermal fluctuations can be rewritten as 

(A) ~ / dsAe ^ , (37) 
Jo 

where A is the torque component that arises from the aver- 
age over free-torque motion and q = 2IiE/J with E is the 
total kinetic energy. It is convenient to define 

A{s,J) = Ae~''^\ (38) 

which represents the torque resulting only from the average 
over free-torque motion as a function of J and the density 
of states in phase space s (see Appendix D). 



In Fig.[22]we show torque components defined by equa- 
tion (|38|) for spheroidal and ellipsoidal AMOs with J = lOJth 
and J = 0.5 Jth, corresponding to cases in which the role of 
thermal fluctuations is marginal and important, respectively. 

It can clearly be seen that for J = lOJth (i.e. J ^ 
Jth), the obtained torques are nearly similar for irregular and 
axisymmetric grains. Yet the torques drop very rapidly to 
zero as s (note that s = sin 7 for the spheroid) increases (see 
dashed lines in Fig. [22]). The former stems from the fact that 
for the suprathermal rotation, a good coupling between the 
maximal inertia axis and the angular velocity is achieved. 
Therefore, there is no difference between the free-motion of 
an irregular and axisymmetric grains. 

However, the averaged torques become very different as 
J decreases. In fact, Fig.[22](see the curves with J — 0.5Jth) 
reveal that for axisymmetric grains, RATs decrease regularly 
with respect to s. Whereas, RATs for irregular grains drop 
rapidly, change the sign and rise again with s increasing. 

The effect of such a variation on the average of RATs 
over thermal fluctuations (see equation I37p is evident. In 
Fig. [23] we show the corresponding torque components (F)^ 
and {//)</) obtained by averaging F and H using equation 
(|37|) . as functions of J for several angles ^ between J and 
k for ■(/; = 0°. It can be seen that for J Jth, the aver- 
aged torques of irregular grains are similar with those of ax- 
isymmetric grains, and they become saturated in both cases. 
However, their averaged torques become very different when 
J decreases. For instance, the averaged torques for the for- 
mer case change their sign at J ~ Jth, while the torques for 
the later case do not. We note that for the angle cos ^ < 0, 
the torque (H)^ is negative (i.e. torques tend to drive grains 
to zero-J attractor point) for J ^ Jth, equivalent with the 
absence of thermal fluctuations. As a result, the change of 
sign of {//)</) for J ~ Jth in the presence of thermal fluctua- 
tions reveals that grains can be spun up again for J ~ Jth- 
To some value of J, (H)^ continues to reverse its sign, and 
grains are decelerated. In other words, the irregular grains 
can be aligned at low attractor points with angular momen- 
tum J ~ Jth- 

Consider now a stationary point ^s. As discussed pre- 
viously, the stationary point is either an attractor or a 
repellor point depending on the sign of the spin-up torque 
and the derivative of the alignment torque at that point; 
it is an attractor point if tttt— ^^^S^ < 0. Therefore, the 

as, ' 

change in sign found above for (F)^ and {H),p can give rise 
to the thermal attractor points observed in maps for the 
ellipsoidal AMO and irregular grains . 

In summary, the analytical analysis for the AMO re- 
sults indicate that as thermal fluctuations become more im- 
portant (i.e. J small), RATs as a function of ^ can have 
different forms (sign and magnitude) that increase the an- 
gular momentum of the low attractor points from J — 
to J ~ Jth ■ We also see a radical difference between torques 
averaged over free precession of spheroidal and the free wob- 
bling of ellipsoidal grains. 

6.2 Stability of low and high-J attractor points 

From Section 6.1, we see that the effect of thermal fluctu- 
ations on RATs is to produce attractor points at low J. 
In addition, there are also attractor points at high J to 
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Figure 22. Spin-up and alignment torques as functions of 5(7) 
for two values of angular momentum J = 20Jth, and J = J^h 
corresponding to an angle t; = 160° between angular momentum 
and the magnetic field, and the precession angle <j> = 180°. The 
fast drop of F and H for small J is observed that gives rise to 
the fact that the sign of their averaged value over s is opposite to 
its sign at s = 0. 



Figure 23. Aligning (-F)0 and spinning (-ff)^ torque components 
averaged over the Larmor precession angle and thermal fluctua- 
tions (i.e. over s using equation I37|l . normalized over their ampli- 
tude values for three angles as functions of angular momentum 
for the spheroid and ellipsoid AMOs. Plots show the change of 
sign of both torques when J'^ / J^y^ decreases (i.e., thermal fluctu- 
ations increase) for the ellipsoid AMO, but not seen for spheroid. 
For J ^ Jtii, the torques are constant. 



which thermal fluctuations have marginal influence. It can 
be shown that most of grains in the ensemble tend to get 
in low-J attractor points. However, having low J, they may 
be significantly influenced by randomization processes (e.g. 
collisions by gaseous atoms). 

For high-J attractor points present in the trajectory 
map (see Figs l9l [141 and [TSl ) . the angular momentum is de- 
termined by the radiation energy density Urad (see equation 
[22)) . that is J /Iilut ~ Urad{H)4,- As a result, depending on 
the distance to a given radiation source, the angular mo- 
mentum of the high-J attractor points changes. 

However, for low-J attractor points, because the angu- 
lar momentum depends on thermal fluctuations (see Hoang 
& Lazarian, in preparation), it does not change when radi- 
ation intensity varies. So grains of the same size, near and 
far from the pumping source, have low-J attractor points of 
similar angular momentum. Therefore, due to thermal wob- 
bling (see Lazarian 1994), even near the radiation source, 
the alignment degree is not high for the case of trajectory 
maps without high-J attractor point. Consequently, even in 



regions close to the strong radiation source (e.g., star form- 
ing regions) most of grains may rotate rather slowly. 

Furthermore, the angular momentum of low attractor 
points depends on the grain size 0^/ / because internal relax- 
ations are a function of Oe// (scales as see also Table 
1). So we expect that for large grains that have weak inter- 
nal fluctuations, the value of angular momentum at low-J 
attractor points is very close to zero as a result of the de- 
celeration action of RATs. In contrast, small grains having 
strong internal relaxations, can have low-J attractor points 
of higher angular momentum. 



7 CROSSOVER DYNAMICS 
T.l Description 

The randomization of grains during crossovers has been 
studied by Spitzer &McGlynn (1979) and Lazarian & Draine 
(1997). In Paper I we studied the crossover for a spheroidal 
AMO. It is shown that grains move to the attractor point 
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cos 5 = — 1, J = while undergoing multiple crossovers. In 
this section, we first consider large grains, so that tc < tsar 
and thermal flipping is not important for both axisymmetric 
and irregular grains. Therefore, we can solve the equations 
of motion for both J, ^ and the component J|j of J along ai. 
Then we incorporate the thermal fluctuations and thermal 
flipping into our treatment. 

In Paper I, to simplify our treatment, we disregarded 
all internal relaxation processes during the crossover. This is 
justifiable as for a^f f 3> Oc the crossover occurs on the time 
scale shorter than the internal (e.g.Barnett and nuclear) re- 
laxation time. However, for typical interstellar grains with 
a^ff < He, thermal flipping occurs very fast, and thus it 
must be taken into account. Because the mirror is assumed 
to be weightless, for the spheroidal AMO, the dynamics of 
free rotation coincides with that of a spheroid. As a result, 
for a given J, ai precesses around J with a constant angle 7. 
The state of the grain is completely determined by describ- 
ing J in the lab and body systems. Equations of motion for 
this case are (see Paper I) 



dJ 

'dt 

Jd cos 7 
di 



J 

tgas 



dJ 



dJn 



, cos 7 + , 
dt ' dt 



(39) 



(40) 



where T is the RAT, and = J cos 7 is the component of 
angular momentum along the maximal inertia axis. 

Assuming that both precession time-scales of ai around 
J and J around B are much smaller than the crossover time, 
we obtain 



§ = f((F>./. + (i^)./-), 



(41) 
(42) 



J d cos 7 
di 



dJ 



{{Qal} + f+ - {Qal)-f-) - — COS 7 - (1 - COS 

(43) 

"^^ , (Qai) is the torque component 



where M — — — 
along the axis ai , /+ , /_ are probability of flnding grains 
in the positive and negative flipping states, respectively (see 
WD03). 

In equation (|43p . the component Qai which is along 
ai also involves the flipping probability because after each 
flipping, 7 TT — 7 for which Qai changes the sign, i.e., 
Qai{T' — 7) ~ —Qai (7) while equation (gSJ describes the 
variation of 7 only due to RATs. 



7.2 Crossover: no averaging over thermal 
fluctuations 

Let us consider first a grain with aeff 3> fflc for which the 
crossover time tc <^ttf. the grain studied is 0.2 micron, and 
this phrase is not correct. In this case, as the thermal hip- 
ping is negligible, the motion of grains comprises the align- 
ment of J with respect to B and the motion of ai about 
J. We solve the equations of motion with the initial condi- 
tion Jo = 2Jth and 70 — n/A. Clearly, when the alignment 
time corresponding to the former case is shorter than the 
crossover time, then the grains are quickly driven to the at- 
tractor points. In particular, grains that become aligned at 



high J attractor points do not undergo crossovers, but oth- 
ers that are driven to low-J attractor points do. The trivial 
consequence of this is that, as 7 changes, it can alter {F}^ 
and (H)^ because they depend on O, which is a function of 
7 and ^ (see Appendix A). 

7.2.1 Only low-J attractor points 

If we consider the crossover of a grains (e.g. the spheroidal 
AMO) for the radiation direction ip — 30° for which only 
low-J attractor points are present in the phase trajectory 
map. The grain dynamics of this grain is similar to the case 
studied in Paper I. 

7.2.2 High- J and low-J attractor points coexist 

Now let us consider a situation when both high-J and low-J 
attractor points are present in the phase trajectory map (e.g. 
tp = 0° for shape 1). As discussed earlier, irregular grains 
undergo irregular wobbling (cf. thermal wobbling) which in- 
duces averaging RATs over three Euler angles. Therefore, 
the obtained RATs are functions of s (see Fig|22p. For sim- 
plicity, we assume that the dependence on 7 of averaged 
RATs for irregular grains is similar to that of axisymmetric 
grains, and thus we can solve equations of motion (|4ip - (|43| l 
for phase trajectories. In fact, the above assumption is fea- 
sible when the slightly irregular grains with I2 ~ /s. 

The results for irregular grains are shown in trajectory 
maps (see Fig. l24p . One result that can be seen directly from 
Fig. [24] is that grains do not reach the low attractor point 
with J = as seen in the phase maps in the absence of 
thermal fiuctuations. Instead, they first tend to reach J < 
Jth, and then reverse their direction to head for the high- 
J attractor point. This stems from the property of RATs. 
Indeed, the torques are functions of O, angle between ai and 
7)J,k. Therefore, we can qualitatively explain the effect based 
on the relationship of ^, O and 7. Firstly, for grains with 
initial angles ^ < 7r/2,7 — 0, the spinning up torques H < 
and F > (see Fig [23]). Therefore, the angular momentum 
of grains decreases while their alignment angle ^ increases. 
In a short time, 7 increases, which gives rise to the change in 
sign of F and H (see also il6.ip . Thus, grains start to reverse 
their orientation (^ decreases) and they are spun up to the 
high-J attractor points. Finally, a perfect alignment at the 
high attractor point is established. 

We see here the difference in the crossover for irregular 
grains and the spheroidal AMO. For example, the difference 
is evident from the different variation of RATs averaged over 
free-torque motion for the spheroidal and irregular grains 
(see Fig. I22p . There it is shown that RATs for irregular 
grains {ellipsoidal AMO) change the sign at small devia- 
tion angle 7 compared to 7 = 7r/3 for the spheroidal AMO. 
It turns out that when grains are driven to low-J attractor 
points, they can be spun up again as seen in Fig. [24] (cf. Fig. 

Fig. [24] (lower panel) represents the phase map for the 
parallel component Jy . It reveals that Jy evolves tightly with 
J. However, grains undergo multiple crossovers (i.e. Jy — 0) 
before being spun up to high J attractor points. Fig. 1241 also 
indicates that during the crossovers, the angular momentum 
of grains is about the thermal angular momentum associated 
with the grain temperature. 
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Figure 24. The phase maps for the angular momentum J ( Upper 
panel) and the component of J along the maximal inertia axis 
ai {Lower panel). The lower panel shows that grains experience 
multiple crossovers. 



One appealing feature seen in the phase map is that 
without thermal fluctuations, grains are not trapped at low- 
J attractor points. They can be spun up after crossovers. 
However, as thermal fluctuations are considered, the results 
in §4 and 5 show that grains are still trapped at low-J at- 
tractor points. 



7.3 Crossover: Averaging over thermal 
fluctuations 

For grains of effective size comparable to the critical size, 
the thermal fluctuations are important because of the de- 
pendence of the Barnett relaxation time on the grain size 
(see Table [!}. In this case, we can consider the grain dy- 
namics in two different regimes. First, assuming that for J 
greater than Jth, the crossover time is shorter than the Bar- 
nett relaxation time, then the crossover can be treated as 
in §7.2. As a result, grains undergo multiple crossovers (see 
Fig.[25]for the phase trajectories with J > Jth)- As J ^ Jth, 
thermal fluctuations occur faster than the crossover time, 
therefore, RATs must be averaged over thermal fluctuations. 
The dynamics of grains in this stage is the same as studied in 
§4 and 5 (i.e. grains are trapped at low-J attractor points). 




Figure 25. The phase map of grains into different regimes corre- 
sponding to the cases with and without thermal fluctuations. For 
J > Jth, this figure represents phase trajectories without averag- 
ing over thermal fluctuations; grains undergo multiple crossovers. 
For J ^ Jth, this figure shows the regime in which RATs are aver- 
aged over thermal fiuctuations; it is shown that grains get aligned 
at the attractor point A, which is similar to what obtained in §6. 



The critical size Oc depends on the radiation energy in 
the anisotropic flux as u^^^ . Thus the alignment near stars 
and in the diffuse ISM is perhaps different. 



8 FAST ALIGNMENT 

In Paper I we showed that the grains can be aligned over 
a time scale much smaller than the gas damping time. Now 
let us consider this problem when thermal wobbling and 
flipping are important. 

Fig. [21] shows the trajectory map constructed from an 
ellipsoidal AMO with QT^'/QS""' = 0.78 and the light 
direction i/) = 70°. The arrow represents a time interval 
At = lOtphot where tphot deflnes the time-scale over which 
RATs decelerate grains from J = Jth to J = in the absence 
of gas damping. It can be seen that the angular momentum 
of the low attractor point is the same as the case in which 
the gas damping is included (cf. Fig. I14[). We observe that a 
significant number of grains are aligned on low-J attractor 
points A, B, C and D over tan ~ LOtphot (see Fig. I26|l . 

As discussed in Paper I, this type of fast alignment can 
occur in a diffuse medium with high radiation intensity, such 
as in the supernovae vicinity or comet wakes. This stems 
from the fact that strong radiation can drive grains faster to 
the low- J attractor points where the thermal fluctuations 
act to maintain the stability of grain alignment. 



9 INFLUENCE OF RANDOM 

BOMBARDMENT BY ATOMIC GAS 

In Paper I we showed that in most cases RATs align grains 
with respect to magnetic flelds while decreasing the grains's 
angular momentum to J = 0. The results above indicate 
that thermal fluctuations can change a zero-J attractor 
point to a J ~ Jth attractor point. In addition, atomic col- 
lisions can affect the alignment established by RATs in the 
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Figure 26. Similar to Fig. 1141 but the gas damping is neglected. 
The arrow represents a time interval 10tpf,ot. 



presence of thermal fluctuations. In this section, we first 
briefly discuss the method of studying gas bombardment 
based on the Langevin equation. Then, we show the influ- 
ence of this process on RAT alignment. 

9.1 Method 

The effects of collisions in the framework of paramagnetic 
alignment have been studied by many authors (Jones & 
Spitzer 1967; PurceU 1971; PurceU & Spitzer 1971; Lazar- 
ian 1997) using the Fokker-Planck equations. However, the 
Langevin approach was used to study this problem numeri- 
cally in Roberge, Degraff & Flatherty (1993) and Roberge & 
Lazarian (1999). The afore-cited papers used the equivalence 
of the Fokker-Planck and the Langevin equations to simulate 
the evolution of angular momentum of grains in a gas co- 
ordinate system. They derived explicit diffusion coefficients 
for random torques acting on spheroidal grains. According 
to the above works, an increment of angular momentum re- 
sulting from gas-grain collisions within an infinitesimal time 
interval dt is (see Roberge et al. 1993) 



dJi = Ai {t)dt + Bij (J, t)dwj , i = x,y, 



(44) 



where dtoj are Wiener coefficients, Ai,Bij are diffusion co- 



efficients, given by 



Ai = (AJi), i = x,y,z, 
(BE''),, = {AJ,AJj),i,j ^x,y,z, 



(45) 
(46) 

where B^ is the transposal matrix of B. Similar to Roberge 
et al. (1993), diffusion coefficients are first calculated in the 
grain body frame, and then transformed to the lab system. 
The diffusion coefficients are averaged over the precession 
of ai around J and over the Larmor precession angle of J 
about B, given by 

(AJ), = (47) 

where tgas is the gaseous damping time. For a spheroidal 
grain, the gaseous damping time is 

O Tb 



igas — 



40r nmvthb^T\\ (em)' 



(48) 



where /^^ is the moment of inertia of the grain along zz axis 
(i.e., along the maximal inertia axis ai), vth ~ ^2mkTgas 
is the thermal velocity of atom, and bm is the semi-axis of 
the grain. r||,rx are factors characterizing the geometry of 
grain given by 



Til = —3 + 4(1 - eS^)fl(e™) - e^^[l - {1 - el 



Fx = 32 - + (1 - e^)5(em) + (1 - 2e„2) 



(l + e-^)2[l-(l-e^)5(e„0]]]. 



(49) 
(50) 



g{em) is related to the eccentricity of the grain through the 
expression: 



, N 1 . ,1 + er, 



(51) 



where Cm ~ \/f — {a/b)'^. Diffusion coefficients are diagonal 
and given by the following expressions in the lab coordi- 
nate system in which the z axis is along the magnetic field 
(Roberge et al. 1993) 



{{AJ^f) = ^nm^btvUl + ^) 

X [(l-Hcos^C)rx -^sin'CTii], 
{{AJy)') = ((AJ.)^), 

((AJ.)^) = ^nmX«L(l+ g) 
X [sin^^Fx -f cos^Cry]. 



(52) 
(53) 



(54) 



Note that the above diffusion coefficients are derived by 
assuming perfect internal alignment of ai with J and for 
spheroidal grains. However, for the sake of simplicity, we can 
adopt these diffusion coefficients for studying the influence 
of gas bombardment on the alignment of irregular grains. 



9.2 Results 

First, we solve the Langevin equation (i.e. equation 144 |l 
for grains subjected to random torques as a result of gas 
collisions, assuming that the diffusion coefficients remain 
the same for irregular grains. We use the initial condition 
J ~ Jth,gas ~ \/2IikTgaa, £, is generated from a uniform 
distribution in the range to tt and (p is a free parameter. 
Wiener coefficients dcu in equation (I44p are generated from a 
Gaussian distribution function at each time-step. Then the 
resulting solution , Jy and Jz are taken as input param- 
eters for solving the equations of motion of grains driven 
by RATs in the spherical coordinate system described by 
J, ^, (see equations (15)-(17)) to obtain new values of J 
and This process is performed over A'^ = 10^ time-step 

At = m^Hn 



Other physical parameters for the ISM are 
taken from Table [T] The alignment angle ^ is used in av- 
eraging over the total time to obtain the degree of external 
alignment, and J is used to calculate the internal alignment, 
assuming that thermal fluctuations follow a Gaussian distri- 
bution. 

For the AMO, we consider our default model of q = 45° 



Here we average over the Larmor precession angle 
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(see Paper I). For DDSCAT, we study the entire spectrum 
of the ISM for grain shapes 1 and 3 with size tie// = 0.2^m. 
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9.2.1 Influence of randomization to the phase trajectory 
map in the presence of high— , J attractor point 

We can see that colhsions have an interesting effect on grain 
dynamics when the high attractor point is present. Random 
collisions are very efficient for low- J. So, the motion of grains 
is disturbed by random collisions when they approach the 
low-J attractor point. After a time interval, grains enter 
the region of cos ^ > (i.e., positive spinning radiative 
torque) for which RATs can spin up to the high J attractor 
point. Therefore, for this case, random collisions increase the 
degree of alignment. The percentage of grains present in the 
vicinity of the high- J attractor points as a function of time 
is shown in the upper panel in Fig. 1271 The respective degree 
of alignment is shown in the lower panel. 

From the upper panel (Fig. [27]), we find that, during the 
time interval t = tgas, only about 10% of grains are present 
at the high — J point. Then it increases with time and attain 
the saturated value of 75% after t = AQtgas ■ 

Following the lower panel in Fig. 1271 we see the rise of 
the degree of alignment with time. It get to a significant 
value after t — lOtgaa, and the alignment is nearly perfect 
with R = 0.8 after t = SOtgas. 

Due to the stochastic properties of gas bombardment, 
the degree of alignment depends on the angular momentum 
of the high attractor point Jhigh{'4')- A detailed study of the 
degree of alignment with JhighO') in Hoang & Lazarian ( ; 
in preparation) shows that the alignment at high attractor 
points is nearly perfect if Jhigh{ip) ^ ■iJth,gas where Jth,gas 
is the angular momentum corresponding to the temperature 
of the ambient gas. 
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Figure 27. For the ISRF, shape 3, and ip = 70°:Upper panel 
shows the variation of the percent of grains present in the vicinity 
of the high J attractor point A, B in time, and lower panel shows 
the Rayleigh reduction factor for this case. 



9.2.2 Influence of randomization to the phase trajectory 
map in absence of high-J attractor point 

The degree of alignment for the spheroidal AMO and irregu- 
lar grains in the absence of high-J attractor points is shown 
in Fig. 1281 For these cases, random collisions act to decrease 
the alignment. This arises from the fact that random colli- 
sions remove grains from the low-J attractor points. We have 
found that if the angular momentum of the high-J station- 
ary point, Jhighiip), larger than Jth,gas, RATs move grains 
to the vicinity of the high-J stationary points, and deceler- 
ate them again. It is seen that although R is decreased by 
gas bombardment, it is still not negligible (e.g., about 0.1 
and 0.2 for AMO and irregular grains; see Fig. I28|) . 

For Jhighi'^P) < Jth,gas, grains are in a fully thermal 
regime. Therefore, the phase map of grains is mostly ran- 
dom. From the lower panel in Fig. 1281 it follows that the 
degree of alignment is marginal. However, a more elaborate 
treatment in Hoang & Lazarian (in preparatiom) does not 
use equation (2) and demonstrate higher degree of align- 
ment. 



10 INFLUENCE OF H2 FORMATION 

TORQUES ON GRAIN ALIGNMENT 

When a hydrogen atom sticks to a grain, it starts its difi^u- 
sion over the grain surface. In general, the grain surface is 
never uniform, and there are always certain special catalytic 
sites where hydrogen atoms can be trapped (see Purcell 
1979; Lazarian 1995). A wandering H atom on the grain sur- 
face may encounter another atom dwelling at the catalytic 
site and a reaction takes place, producing a II2 molecule. 
The ejected H2 molecule acts as miniature rocket thrusters. 
Averaging over the grain surface (e.g. a brick surface), H2 
rockets produce a net angular torque that is parallel to the 
maximal inertia axis (Purcell 1979). However, resurfacing 
or poisoning by accretion of heavy elements can destroy the 
catalytic sites and create new ones. As a result, H2 torques 
change both magnitude and direction over a definite time 
scale tL, so-called resurfacing time. 

10.1 Method 

Purcell (1979) used the Monte Carlo method to simulate 
the variation of the torque as a result of grain resurfacing. 
Roberge & Ford (2000) used the equivalence of the Langevin 
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Figure 28. Figs show dynamics of external, internal degree of 
alignment and Rayleigh reduction factor in time corresponding to 
RATs from the AMO {upper panel) and DDSCAT {lower panel) 
when the phase trajectory maps do not have high-J attractor 
points. 



and the Fokker-Plank equations to simulate the H2 torque. 
They model Lz as a Gaussian process that is averaged to 
zero in time, and the correlation function that is exponential: 



LUt)LUt - r) = ((AL,)")exp(-r/ti). 



(55) 
(56) 



In equation (|56l) . {(AL^)^) is the magnitude of H2 torque 
given by 

{(ALz)^^' = ^-('^f''MH){2EkT,asf^\%fl, (57) 

where / is the fraction of H atoms absorbed by the grain 
and converted to H2, n(H) is the H density, E is the kinetic 
energy of each departing H2 and I is the side of the individual 
catalytic site (Purcell 1979; Lazarian & Draine 1997). 

In an interval of time dt, the torque along the axis ai , 
, can be simulated by a Gaussian process with the correlation 
time scale tL, given by the Langevin equation (Roberge & 
Lazarian 1999) 

dLl^~L':- + [2^^^^]'^^dw, (58) 
where dw is a random variable, independent of time, sam- 



pled from a Gaussian distribution function, and L'!,{t) is the 
instant torque at time t. This equation allows us to follow 
the evolution of H2 torque in time. 

When the internal relaxations (Purcell 1979; Lazarian 
& Draine 1999b) are taken into account, the angle 7 be- 
tween the maximal inertia axis ai and the angular momen- 
tum fluctuates in time. However, as the fluctuation time- 
scale is much shorter than the alignment time-scale, we do 
not follow the evolution of this angle. Instead, we average 
over this. Following Spitzer & McGlynn (1979), the mean 
H2 torque for a spheroid is given by 



= Lz cos 7 J, 



(59) 



where cos 7 denotes the average of cos 7 over thermal fluc- 
tuations as defined by equation p2p . which defines the value 
of the projection of onto the angular momentum axis 
J\\Z. For irregular grains, cos 7 in equation ((59} is replaced 
by the thermal average, C, of C which is given by 



h - /3'7 sl/2 



and 



\l/2 



h 



2F{Ti/2,k- 



-, for q > 



h 



(60) 



(61) 



where and F are given in Appendix C (see WD03). 
Because of the grain flipping, equation (|59|l becomes 



Th2 = L\ cos 7(/4 



(62) 



where /+ and /_ are probability of flnding the grain in pos- 
itive and negative flipping states, respectively. 

The equations of motion (i.e. equations 1 16m7|) of grains 
by RATs and H2 torque read 



di ^ 7^^rad"eff^ 
dt 
dJ 



2J 



dt 



-^lUj.g^d"-lgX{H{^, (j), V, J)}4> - 



to 



H2: 



(63) 
(64) 



where we have averaged RAT components over the Larmor 
precession angle (p. 



10.2 Results 

We solve the equations of motion for grains subjected to 
RATs and H2 torques using the Runge-Kutta method with a 
finite time-step. At — 10~'''tgas and = 10^ time-steps . We 
use the initial condition J — 20 Jth, and ^ is generated from a 
uniform angle distribution in the range from to pi. At each 
time-step, we first solve the Langevin equation (equation 
I58p for components of H2 torques. Then, by subsituting the 
resulting H2 torques into equations I63|l - (|64|l . we solve the 
obtained equations for J and ^. The resulting solutions J 
and ^ are used to construct trajectory maps and calculate 
the degree of alignment. 



10.2.1 Trajectory maps 

Fig. [29] shows the map when the variation of H2 torques is 
accounted for (lower panel), compared with the map driven 
by RATs (upper panel) produced by a radiation field of 
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Figure 29. Upper panel: map trajectory of grain Qe// = 0.2fj,m 
and A = 1.2fim by RATs; lower panel: when H2 torque varies 
due to resurfacing of active sites on grains surface. The results 
show that, during the spin-down, grains are driven by RATs to 
thermal angular momentum. However, because grains flip very 
fast, H2 torques become inefficient to spin up grains at low-J 
attractor points. 



A = 1.2^m for a grain with size (Je// = 0.2^m. We as- 
sume that the resurfacing process on the grain surface oc- 
curs rapidly, so we model the variation of H2 torques with 
the correlation time-scale tL = O.ltgas - In the absence of H2 
formation torques, RATs drive grains to attractor points. 
A significant fraction of grains are driven to low attractor 
points marked by C, D, at J/Iiujt ~ 1, and some grains 
are aligned at attractor points A, B corresponding to high 
angular momentum (see the lower panel in Fig. I29|l . 

Consider now a case of when H2 torque amplitude is 
larger than that of RATs (e.g. consider the radiation direc- 
tion if) — 70°). It is obvious that RATs still dominate the 
Purcell torque for this case, and drive a significant number 
of grains to the low attractor points C and D. In fact, at the 
initial phase of J ^ Jth, grains are spun up by the Purcell 
torque; however, during the spin-down, RATs lead grains to 
low-J attractor points. Because of grains having low angu- 
lar momentum at the low-J attractor points, grains fiip fast 
and the Purcell torques are averaged to zero (see the lower 
panel in Fig. I29|l similar to what is described in Lazarian 



& Draine (1999a). Therefore, the trajectory map with two 
attractor points C, D of low angular momentum is deter- 
mined by RATs as in the case without H2 torques (see the 
upper panel in Fig. I29p . However, for the attractor points 
A, B with high angular momenta, the thermal flipping does 
not occurs because of J S> Jth- In other words, during spin- 
down period, grains are driven by H2 torques to reach the 
low angular velocity regime. As J decreases, RATs become 
dominant and drive grains to attractor points C, D, as shown 
in the lower panel of Fig. (55] However, as J increases, H2 for- 
mation torques become dominant and do not allow grains 
to have attractor points. Therefore, attractor points A, B 
become unstable, and grains are spin up and down by H2 
torques and RATs in the range between X and Y (see the 
lower panel in Fig. [29]), corresponding to J /Jth = 5 to 40. 
The range of such a variation depends on the relative am- 
plitude of J/2. 

When the correlation time-scale II = Wtgas, our results 
show that trajectories of grains ending at high- J attractor 
points C, D are marginally affected by H2 torques, because 
RATs eventually drive grains to these attractor points at 
which the fast thermal flipping destroys completely the in- 
fluence of H2 torques. In contrast, grains that are driven to X 
and Y experience the same effect as in the case t_t ~ O.ltgas- 

For II —> 00, H2 torques act to spin grains up. There- 
fore, the low-J attractor points become high -J attractor 
points, provided that H2 torques dominate RATs. 

For larger grains, thermal flipping is less important (see 
Lazarian & Draine 1999a) , so H2 torques are important even 
during the slow rotation period, influencing significantly the 
radiative alignment. We provide the treatment of this case 
elsewhere. 



10.2.2 Dynamics of degree of grain alignment 

Fig.[30]shows the variation of the degree of alignment (i?) in 
time for the cases in which only RATs are considered (upper 
panel) and RATs plus H2 torques with the correlation time- 
scale tL ~ O.ltgas are taken into account (lower panel). 

It is seen from Fig. |3D] that the degree of alignment (R) 
is the same for the case ?/> = 70'^ because of the dominance 
of RATs in driving grains to attractor points. At ip = 30", 
the alignment degree exhibits small variations. In addition, 
for both directions, (R) does not increases smoothly to a 
stable value as in the case of alignment by only RATs as 
grain dynamics is affected by the stochastic variation of H2 
torques. When grains are driven to the attractor points with 
thermal angular momentum A and B, the fast thermal hip- 
ping of grain wipes out the effect of H2 torques, which are 
fixed in the grain body coordinate system. Therefore, the 
grain alignment is mostly determined by RATs. 

10.3 Alignment by RATs, H2 formation torques 
and gas bombardment 

In most cases, RATs play the key role in aligning grains 
with a magnetic field, while H2 is only important to spin up 
grains. Moreover, for grains smaller than ~ 10~''cm, thermal 
flipping is very fast, and thus the influence of H2 torques is 
rather marginal. In contrast, for large grains ^ 10~* cm, H2 
torques could be significant for the grain alignment . On the 
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Figure 30. Panels show the variation of degree of alignment (R), 
of J with B in time for the monochromatic light illuminating grain 
in two directions = 30° and 70° . upper panel: only RATs align 
grains and they have stable orientation; lower panel: the same 
as the lower panel, but H2 torques with a correlation time-scale 
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Figure 31. The Rayleigh reduction factor for grain a^ff = 
0.2/nn aligned by RATs, H2 formation torques (upper panel) and 
by RATs, H2 and gas bombardment [lower panel). 



other hand, random coUisions affect the alignment of grains 
mostly when grains rotate at thermal velocity. The upper 
panel in Fig. 1311 shows both internal and external degree of 
alignment as well as the Rayleigh reduction factor for the 
alignment induced by RATs and H2 torques. It is seen that 
R ~ 0.25 for this case. However, when gas bombardment is 
taken into account, we found that R increases to 7i = 0.75 
(see the lower panel in Fig. I3ip . In fact, this increase is 
not unexpected, since for this case, grains are driven by gas 
bombardment to high attractor point, which increases the 
degree of alignment as discussed in §8. 



11 DISCUSSION 

11.1 Central role of the AMO 

Our understanding of the RAT alignment has been improved 
considerably when a simple model of a helical grain (i.e. the 
AMO), was introduced in Paper I. With the assumption of 
perfect internal alignment, the AMO reproduced the align- 



ment effects similar to irregular grains studied by DDSCAT 
(see Paper I). For instance, the AMO shows that for a given 
shape, size and radiation field, some grains can be spun 
up to suprathermal rotation and aligned at high-J attrac- 
tor points, while most of grains are driven to low-J attrac- 
tor points. Moreover, the AMO demonstrated that whether 
RATs can spin-up grains to high J attractor points depends 
on the ratio of the torque components Qei and Qe2- 

The correspondence of the AMO to DDSCAT calcula- 
tions of RATs for irregular grains in Paper I encouraged 
us to use the AMO for the case when thermal wobbling is 
taken into account. In our study, we used the AMO assuming 
that its inertial properties are defined by a triaxial ellipsoid 
rather than a spheroid. This was very easy to accomplish, 
as, within the AMO, torques acting upon either spheroidal, 
or ellipsoidal body induce only precession, while the actually 
important torques (i.e. Qei and Qe2) arise from the mirror. 
We found a number of differences in the AMO dynamics for 
the two cases. We treat the AMO with ellipsoid of inertia as 
our generic case. 
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11.2 Treatment of free wobbling 

In Paper I we showed that Qes in the AMO can be neghgible 
when considering the alignment and spin-up torques. We 
have found again that the component Qes for the AMO also 
plays a marginal role in alignment, as well as, spinning up 
grains if wo followed the algorithm in WD03 and averaged 
RATs over a sufficient time-steps (e.g. about A'' = 10^; see 
§3.3). Using our new algorithm of averaging in which we 
solve the dilfcrcntial equation for the rotation angle instead 
of assuming that this angle ranges from to 2tv, wc obtained 
similar results when the averaging is perfomed over a much 
longer time (e.g. P = W^Pt). It contributes mainly to the 
precession about the radiation beam axis k (see Fig. 1). As 
a result, in most practical situations, it is sufficient to deal 
with only two RAT components Qei and Qe2. 

We studied the alignment for two cases (1) when ther- 
mal relaxation time-scale is much shorter than the crossover 
time (i.e., averaging RATs over thermal fluctuations) and (2) 
when the thermal relaxation time scale is much longer than 
the crossover time. For the former case, grains are directly 
aligned at low-J attractor points. In contrast, the grains un- 
dergo multiple crossovers in the second case. 

The AMO enables us to average RATs analytically over 
thermal fluctuations assuming a Gaussian distribution. The 
obtained expressions allowed us to gain important insights 
into the role of thermal wobbling. 

11.3 Direction of grain alignment 

It has been shown in Paper 1 that the alignment by RATs 
occurs mainly with the longer axis perpendicular to the mag- 
netic field. However, there was still a narrow range of radi- 
ation directions for which grains were aligned with longer 
axes parallel to magnetic field (i.e. "wrong" alignment). We 
have shown that the "wrong" alignment can no longer exist 
because of the fast wobbling of the grains. Indeed, in Paper 
I we found that "wrong" alignment takes place for a low-J 
attractor point for a narrow range of angles around the an- 
gle between the light direction and magnetic field ip = ■k/2 . 
This range is narrower than the range of grain wobbling at 
the low-J attractor point. The disappearance of the "wrong" 
alignment was, in fact, predicted in Paper I. A direct impli- 
cation this effect is that allows a more reliable interpretation 
of the observation data: grains are aligned with long axes 
perpendicular to magnetic field when they are aligned by 
RATs. 

11.4 Degree of grain alignment and critical size 

The degree of external alignment of angular momentum with 
the magnetic field can achieve unity. However, because of the 
weak coupling of the maximal inertia axis of grain with the 
angular momentum at low-J attractor points, the Raylcigh 
reduction factor R is lower than unity. A new effect that we 
discuss in this paper is the transition of grains from low-J 
to high-./ attractor points, in situations when such attractor 
points coexist. This increases the degree of alignment. Inter- 
estingly enough, the transition which makes grains better 
aligned is induced by random bombardment (see §9). 

The degree of alignment depends on the angular mo- 
mentum Jhigh of the high attractor points. The value of 



Jhigh is a function of RATs that increases with grain size. 
Because of the gas bombardment, only grains aligned with 
Jhigh > Jth cau maintain the stable alignment with the 
magnetic field. The minimal size corresponding to Jhigh,min 
is named the critical size. For the irregular shape 1 and 
the ISRF, by studying the alignment for grains with size 
spanning from 0.025/xm to 0.2/iTn, we found that the crit- 
ical size of aligned grains is about 0.05^m, assuming that 
Jhigh,min = 3Jth. Apparently, critical size depends on radi- 
ation field. Therefore, for attenuated field, e.g. dark molec- 
ular clouds, the critical size is shifted to larger value (i.e., 
only large grains can be aligned by RATs). 



11.5 Role of paramagnetic dissipation 

Paramagnetic dissipation does not play any role in the RAT 

alignment. Although the predictions of RAT alignment co- 
incides with the predictions by the Davis-Greenstein mech- 
anism (i.e. the alignment with long axes perpendicular to 
the magnetic field) , RAT alignment occurs over shorter time 
scale compared to paramagnetic relaxation time too, pro- 
vided that grains are paramagnetic. 

As a result, unless grains are superparamagnetic (see 
Jones & Spitzer 1967; Mathis 1986), paramagnetic relax- 
ation is irrelevant. A common fallacy entranced in the liter- 
ature is that RAT alignment is a sort of paramagnetic align- 
ment of suprathermally rotating grains, i.e. kind of align- 
ment of fast rotating suggested by Purcell (1979), but with 
fast alignment produced by RATs. This way of thinking of 
RAT alignment is erroneous, as, unlike the Purcell (1979) 
torques, RATs induce their own alignment, which is faster 
than the paramagnetic one. 



11.6 Effects of radiation field intensity 

When grains are aligned at high- J attractor points, the vari- 
ations of radiation intensity induce the variation in J value. 
However, a significant fraction of grains that align at low 
-J attractor points at which J does not vary with radia- 
tion intensity. Therefore, after a certain threshold at which 
grains are aligned suprathermally (i.e. with u) ^ (^thermal) 
at high-J attractor point the increase of radiation intensity 
does not increase the degree of alignment for grains aligned 
at such points. 

When grains are aligned at low-J attractor points, a 
weaker alignment is expected, because of the poor coupling 
of the maximal inertia axis with angular momentum for low 
J (see Lazarian 1994; Lazarian & Roberge 1997), irrespec- 
tive of the intensity of radiation field. The increase of the 
radiation intensity does not affect the position of lower- J 
attractor points. 



11.7 Role of gas bombardment and H2 torques 

This is the first study to take into account the effects of 
gaseous bombardment and H2 formation within the RAT 
alignment mechanism. Random gas bombardment itself has 
always been considered as the cause of grain randomization. 
However, in the framework of the RAT alignment, its effect 
can be different. 
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As discussed in the paper, RATs align grains with re- 
spect to the magnetic field while driving grains to high- J 
and low- J attractor points. Naturally, the later case is more 
sensitive to gas bombardment. This random process can in- 
crease the alignment degree by driving grains from the low- J 
attractor points to high J attractor points. Such a process 
gives rise to the time-dependent alignment. The time over 
which the alignment can become stable is about lO^fgas- For 
high- J attractor points, the alignment is marginally affected 
by the bombardment. 

H2 torques were believed to play a major role on the 
grain alignment. However, we found in the present study 
that H2 torques are indeed important for high attractor 
points for which they can spin grains up. For low attrac- 
tor points with low angular momentum, the fast flipping 
of grain axis as a result of thermal fluctuations averages 
out any torques which are fixed in the grain body, e.g. H2 
torques. 



11.8 Time-scales of alignment 

If paramagnetic dissipation is neglected, one may guess that 
the only characteristic time that can determine RAT align- 
ment could be tgaa (see DW97). As pointed out in Paper 
I, strong radiation field can provide grain alignment over 
a time-scale shorter than the gas damping time (i.e. fast 
alignment). The characteristic time of this alignment can be 
~ SOtphot, where tphot is the time over which the amplitude 
value of radiative torques can deposit a grain with the angu- 
lar momentum of the order of its initial angular momentum. 
The large coefficient in front of tphot reflects relative ineffi- 
ciency of radiative torques in the vicinity of low- J attractor 
points. The supernova shell or star forming regions are fa- 
vorable medium for such a fast alignment. 

Our present study shows that another time scale is ap- 
propriate, when the phase trajectory map contains both low- 
J and high-J attractor points. This time is about ~ SOtgas 
and it corresponds to the time scale over which the gaseous 
bombardment transfers grains from low- J to high-J attrac- 
tor points. 



11.9 Role of initial J 

The role of initial J for alignment has been studied in Pa- 
per I, where we found that it determines the time of "fast" 
alignment tphot- In addition, in the paper above, we have 
seen new low- J attractor points present in the phase maps 
when thermal fiuctuations are accounted for. However, if 
grains start from small J, the grain trajectories may not en- 
ter the parameter space, that would induce directing grains 
to high-J attractor points. Instead, the trajectories may, in 
the absence of gaseous bombardment, end up at the low-J 
attractor points. 

However, in practical terms, this particular difference 
is not crucial for the overall dynamics. In the presence of 
gaseous bombardment, the grains will be moved to the high 
attractor points, if such points exist. 



11.10 Grain dynamics at J ~ Jth 

In Paper I, with the assumption of simplified dynamics, i.e., 
assuming that J is always aligned with the maximal inertia 
axis ai , we showed that AMO induces a good alignment with 
respect to radiation direction or magnetic field. The assump- 
tion of J and ai coupling is correct only for J 2> Jth- In this 
paper, we use the analytical formulae for RATs and get the 
analytical expressions for torques components when thermal 
fiuctuations are accounted for, and therefore the angle be- 
tween J and ai fiuctuates. Results show that torque compo- 
nents decrease substantially as thermal fiuctuations become 
stronger (see Figs [7]), but the dynamics of grains does not 
change radically. This interesting finding testifies that ther- 
mal fiuctuations as a result of internal relaxations do not 
play the key role in creating the low-J attractor points as it 
was believed earlier, but they "lift up" the earlier existing 
attractor points. 

In order to see whether the irregularity of grain shape 
is important, we modify slightly AMO by replacing the 
spheroidal body by the ellipsoidal one. As earlier discussed, 
the eUipsoidal body does not affect on Qei,Qe2 for AMO, 
but changes the dynamics of fiuctuations and therefore the 
averaging. As a result, the torques resulting from averag- 
ing over thermal wobbling are different, which induces the 
difference in the dynamics for irregular grains. Thus, the ir- 
regularity of shape and the presence of thermal fiuctuations 
act together when grain angular momentum is comparable 
with Jth- 

11.11 Crossovers 

We studied the crossover by following the variations of an- 
gular momentum and the component parallel to the axis of 
major inertia in the case when thermal fiuctuations are neg- 
ligible. Our detailed treatment of crossovers is in agreement 
with the earlier studies by Spitzer & McGlynn (1979) and 
Lazarian & Draine (1997), but is different from the treat- 
ment of crossovers in DW97. 

Whereas in Paper I we followed the regular phase trajec- 
tories of grains approaching to the crossover, in the present 
study we take into account thermal wobbling that becomes 
important as the grain angular momentum becomes compa- 
rable with the thermal one. Unlike the approach in WD03, 
which is applicable only to small grains, our treatment is 
also applicable to large grains for which initially the time 
of thermal wobbling (which is of the order of the internal 
relaxation time) is longer than that of crossovers. We find 
that both large and small grains eventually get to a similar 
state, at which the amplitude of the angular momentum gets 
of the order of its thermal value. This corresponds to low-J 
attractor point. 

Our work confirms the claim in Paper I that the low-J 
attractor points are not metamorphosis of crossover points 
after which grains spin up again as in DW97 (see also 
WD03), but a modification of zero- J attractor points, that 
exists in the absence of thermal fiuctuations. 

11.12 Helicity of grains 

Helicity of grains as a factor of RAT processes was first men- 
tioned in Dolginov & Mytrophanov (1976). However, there it 
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was not identified what can make grains helical. Moreover, 
it was claimed that some regions of space can have non- 
helical grains, which implied that helicity of grains is some- 
thing that requires rather special conditions to emerg 
Later the RAT action was demonstrated for rather arbitrar- 
ily shaped grains in DW96, but this and subsequent studies, 
namely, DW97, WD03, did not use the concept of helicity. 

Grain helicity, as an essential requirement for alignment 
was identified in Paper I, where by comparing results of nu- 
merical calculations for irregular grains with AMO, it was es- 
tablished that grains can be classified as either of left or right 
helicity. Moreover, it was shown that irregular grains demon- 
strate helicity and, therefore, alignment, not only subject to 
radiation, but also to gaseous flows (see more in Lazarian & 
Hoang 2007b). 

In Paper I, an irregular grain rotating about its axis 
corresponding to the maximal inertia axis was identified as a 
helical grain. Our present study demonstrates the advantage 
of good internal relaxation of the maximal inertia axis and 
angular momentum. Indeed, we show that torques acting on 
a grain decrease when the grain starts wobbling. 

11.13 Irregularity of shape and inertia 

Irregular grains, in general, have an irregular shape and have 
to be characterized by a triaxial ellipsoid of inertia, rather 
than a spheroid. The irregularity of grain shape gives rise to 
grain helicity, while the irregularity in terms of the moment 
of inertia modifies grain dynamics. One of the interesting 
effects that we observed for grains, which can be character- 
ized by a spheroid of inertia, is a transfer of grains from 
low-J attractor point (actually, zero-J attractor point) to 
high-J attractor point in the absence of gaseous bombard- 
ment. This effect disappears when thermal fluctuations are 
accounted for. 

11.14 Comparison with earlier studies 

The irregular torques arising from emission and absorption 
photons randomly depositing angular momentum to grains 
were discussed in Harwit (1970). The marginal importance 
of the mechanism for the alignment was shown, however, in 
Purcell & Spitzer (1971). Dolginov (1972) proposed the first 
model of RATs that act on chiral, e.g. quartz, grains. Later, 
Dolginov & Mytrophanov (1976) considered a more generic 
case of RATs, i.e. RATs arising from "twisted", irregular 
grains. The idea of RAT alignment was mostly ignored by 
the community (see Lazarian 1995, as an exception) until nu- 
merical calculations of RATs got available in DW96, DW97. 
From the earlier studies, the most relevant to this work is 
WD03, where effects of thermal fluctuations within the RAT 
alignment were empirically studied. 

This work extends our study in Paper I, which demon- 
strated the ability of the AMO to represent helical grains. 

An additional confusion in the aforementioned study stemmed 
from their claim that the prolate grains should align with long 
axes perpendicular to the alignment of the long axes of oblate 
grains. However, it was shown in Lazarian (1995) that both types 
of grains align the same way, if internal relaxation is accounted 
for. 



The AMO allowed us to address the generic features of the 
RAT alignment. We note that our procedures of torque av- 
eraging and our results for angular momentum for low-J 
attractor points differ from those in WD03. Our interpreta- 
tion of the origin of those points is also different. 

Extending AMO, we have proven that thermal trapping 
of grains at the low J attractor points reported in WD03 
is generic, i.e. does not depend on the particular direction 
between the beam and the magnetic fleld, or on the partic- 
ular choice of the irregular grain or on the particular wave- 
length/radiation spectrum chosen. Moreover, we have shown 
that the we expect the new low J attractor points to appear 
close to cos^ = 1, which makes grains aligned with B at 
these new attractor points. All the results we obtained with 
AMO were also tested with two irregular grain shapes. In 
addition, we accounted for gaseous bombardment and H2 
formation. 



11.15 Accomplishments and limitations of the 
present study 

The study above clarifled a number of outstanding issues 
in RAT alignment theory. It relied on the guidance of the 
analytical studies of RATs and crossovers in Paper I and 
accounted or the proper averaging of RATs arising from 
thermal wobbling. It confirmed the predicted in Paper I 
the gradual change of the position of lower attractor points 
on the trajectory maps of irregular grains from zero to the 
thermal value of angular momentum. In addition, our study 
accounted for the effects of gaseous bombardment and II2 
formation on the RATs. This resulted in a discovery of a new 
important effect, namely, the transfer of grains from lower 
attractor points to high attractor points, when the latter 
points exist. This is a single most important finding of our 
paper. It is also a practically important effect, as the internal 
alignment is perfect for grains at high attractor points and 
reduced for grains at lower attractor points. All in all, the 
present paper substantially extends our analysis in Paper I. 

The limitations of the present study stem from the fact, 
that, while addressing the issue of how the degree of align- 
ment changes in different circumstances, it does not calcu- 
late the degree of alignment exactly. This is done in Hoang 
& Lazarian (in preparation). 



11.16 Towards modeling of grain alignment 

Present day modeling of grain alignment both in molecular 
clouds (CL05; Pelkonen et al. 2007; BetheU et al. 2007) and 
accretion disks (Cho & Lazarian 2007) uses rather heuris- 
tic recipes for determining whether grains are aligned. In 
these studies the amplitude values of the RATs are calcu- 
lated and used to calculate the maximal angular velocities 
achievable for a given damping of grain rotation. Such ve- 
locities parametrize the efficiencies of grain alignment in the 
chosen environments. In the paper above, we have confirmed 
the utility of such a parametrization and improved the crite- 
rion for the alignment to be efficient. However, we still have 
to obtain better measures of the expected alignment. 
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12 SUMMARY 

In the present paper, we studied the role of thermal fluctua- 
tions, thermal flipping, efficiency of H2 torques and influence 
of random collisions on the RAT grain alignment for both 
the AMO and irregular grains. Our main results are sum- 
marized below: 

1. We have used an AMO with inertia defined by a tri- 
axial ellipsoid and averaged numerically RATs over thermal 
fluctuation. For this ellipsoidal AMO, we have found the in- 
crease of angular momentum of low- J attractor points from 
J = to J ~ Jth- We also found the similar effect for ir- 
regular grains in which RATs are calculated using DDSCAT 
and for the entire spectrum of the ISRF. 

2. We have proved that "wrong" alignment correspond- 
ing to low angular momentum for a narrow range of radia- 
tion direction reported iu Paper I is eliminated when thermal 
fluctuations and thermal flipping are considered. 

3. We have studied the crossover for the AMO and ir- 
regular grains and found that grains experience multiple 
crossovers before being trapped at low-J attractor points 
in the presence of thermal fluctuations. 

4. We have found that random collisions of atomic gas 
increase the degree of alignment when grains are aligned by 
RATs with both low-J and high-J attractor points by driv- 
ing grains from low-J to high-J attractor points. When the 
angular momentum of the high-J attractor point Jhigh{'<P) 
greater than 3Jth,gas, a significant degree of alignment can 
be achieved. 

5. We have studied the influence of H2 formation 
torques in the frame of the RAT alignment for different 
resurfacing times, showing that it can, in particular circum- 
stances, enhance the degree of alignment with respect to the 
magnetic fleld. 
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APPENDIX A: RELATION OF ANGLES 9, /?, $ 

AND 5,0,7/) 



cos B = ai.ei. 



$ = 2tan 



if sin 9 — ai.e2 



ai.eg 



(Al) 
(A2) 



P = 2tan ^ 



sin 6 + 32.62 



sin 0(a2.e3 cos $ — 3.2.62 sin $) 



(A3) 

where ai, 32,33 are unit vectors in grain reference, and 
61, 62, 63 are unit vectors in the lab system. Dot products in 
above equations can be obtained from following expressions 

3i.ei = cos tpZB-SLi — sin ijjXB-SLi, (A4) 

3^.62 = COS tpHB-Sii — sin ^ZB-SLi, (A5) 

SLi.e-i =yB-SLi, (A6) 

where 

xs.31 = COS (j> sin ( sin 7 -I- sin (j) cos ( sin 7 

+ sin 5 COS (j) COS 7, (A7) 

zb.3i = cos 5 cos 7 — sin ^ sin ( sin 7, (A8) 

ys.si = cos 5 sin <j) sin sin 7 — cos <j) cos sin 7 

+ sin 5 sin cos 7, (A9) 

xs.32 = cos 5 cos (j>{ cos a cos ^ — sin a sin cos 7) 
— sin (j>{ cos Q sin ^ + sin a cos ^ cos 7) 
-1- sin 5 cos (j) sin a sin 7, (AlO) 

zb.32 = — sin 5( cos a cos ^ — sin a sin ( cos 7) 

+ cos 5 sin Q sin 7, (A-11) 

ys.32 = cos 5 sin <j){ cos a cos — sin a sin ( cos 7) 
-f cos </>( cos a sin ^ + sin a cos ^ cos 7) 
+ sin 5 sin <j!> sin a sin 7, (A12) 

(A13) 

with x_B , ys , z s are unit vectors of magnetic field coordinate 
system in which zb||B. Here Q,7,C are Euler angles shown 
in Fig. 



APPENDIX B: RATS 

Similar to Paper I, in order to make an easy relation of our 
results to those in earlier works, wherever it is possible, we 
preserve notations adopted in DW97. Mean radiative torque 
efficiency over wavelengths, Q(0, /3, "l") is defined as 



— _ / QxUxdX 
J u\dX 



(Bl) 



where u\ is the energy density (see Mathis, Mezger & Pana- 
gia 1983), and Qa is the RAT efficiency corresponding to 
wavelength A. RAT from the anisotropic component of ra- 
diation is defined by 



— 2 T 



-7Q 



(B2) 



where 7 is the anisotropy degree of radiation, a^ff is the 
effective size of the grain (see DW96; Paper I), X,Urad are 
mean wavelength and mean energy density of radiation field, 
which are respectively given by 



A = 



Urad 



J XdX 

J u\XdX 
JXdX ' 



(B3) 
(B4) 



APPENDIX C: AVERAGING OVER 
FREE-TORQUE MOTION 

Consider an ellipsoid with three principal axes 3i , 32 , 33 and 
moment of inertia Ii > I2 > I3, respectively. Dynamics of 
such an ellipsoid is clearly presented in the textbook of Lan- 
dau & Liftshift (1972). Let define a dimensionless quantity 



{h-h)(l-hq/hy 



(CI) 



where 



'iliU is the ratio of total kinetic energy to the 



1 - J 

rotational energy along the maximal inertia axis. 
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For q < /1//2, fc^ < 1, the solution of Euler equations is 
coi=±^{^i-^)'^'dn{r), (C2) 
'''^sn(r), (C3) 



-12 



/3 



/l- 




II - 


h ' 




-1) 


h - 


-h 




-1) 


h - 


-h 



(C4) 



where cn, sn, tan are hyperbolic trigonometric functions, 
and T is given by 



tJ 



{h-h){l-hq/h) 



1/2 



(C5) 



/1/2/3 

The rotation period around the maximal inertia axis is 



Pr ^'iF{-K/2,e), 

where F is the elliptic integral defined by 

F{e,m)=[ de(l-msin^6l)"^''^ 
Jo 

For q > /1//2, angular velocities are given by 

J fh- hq si/2 



ui — 



U12 



cn(T), 



J 12(1 - l3q).i/2 I N 



0^3 =±— ( 



J" ,h[.q - i)ni/2 , , X 



Is ' /l - /3 



where 



(J2-J3)(g-1) 



1/2 



/1/2/3 

Rotation period for this case is given by 
Pr = 4F(7r/2,fc^). 



(C6) 
(C7) 

(C8) 
(C9) 
(CIO) 

(Cll) 
(C12) 



In above equations, ± stand for the positive and negative 
flipping states. 

Euler angles Q, 7, C (see Fig. IA1|) can be deduced from 
angular velocity as foUowings 



sm a sm 7 ; 
cos a sin 7 - 
cos 7 : 



J ' 



(C13) 
(C14) 
(C15) 



WD04 averaged RATs over free torque motion as following. 



^T- Jo 



dr— I A{aX,l). 



(C16) 



Here they assume that the rotation of the grain about the 
angular momentum is much faster than the rotation about 
the maximal inertia axis (i.e. at each a moment t with an- 
gles 7 and a, the angle can vary from to 2-k). Clearly, 
this assumption is not correct because these processes are 
comparable. We have implemented an algorithm to average 
RATs over free torque motion in which the angle ^ is ob- 
tained by solving an Euler differential equation. We found 
that RATs obtained by our method are in good agreement 
with what obtained using equation (C16) when they are av- 
eraged over a time scale greater than IO'^Pt. 



APPENDIX D: AVERAGING OVER THERMAL 
FLUCTUATIONS 

Average of a quantity A over thermal fluctuations is deflned 
by 

_ J^dsAiq,J)eM-q{s)JV2hkU 
^ ' J dseM-qjy^hkTd] ■ ^ ' 

Number of state s in equation ()D1|I is given by 



s = 1 



da 



Jo 

with ai is given by 



/3(/i - hq) + hih - h)cos^a 



[ h{Ii - I2) + hih - h)co8^a 



1/2 



(D2) 



ai = cos 



h{hq~h) 



11/2 



/l(/2-/3) 

for q > h/h, and qi = 7r/2 for q ^ h/h (see WD03) 



(D3) 
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